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Abstract. In this paper, we examine the non-relativistic stationary Schrodinger 
equation from a differential Galois-theoretic perspective. The main algorith- 
mic tools are pullbacks of second order ordinary linear differential operators, 
so as to achieve rational function coefficients ("algebrization"), and Kovacic's 
algorithm for solving the resulting equations. In particular, we use this Ga- 
loisian approach to analyze Darboux transformations, Crum iterations and 
supcrsymmetric quantum mechanics. We obtain the ground states, eigenval- 
ues, eigenfunctions, eigenstatcs and differential Galois groups of a large class 
of Schrodinger equations, e.g. those with exactly solvable and shape invariant 
potentials (the terms are defined within). Finally, we introduce a method for 
determining when exact solvability is possible. 
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1. Introduction 

In classical mechanics, the Hamiltonian corresponding to the energy (kinetic 
plus potential) is given by 

H = + U(~x*), ~f = (pi, . . . ,p n ), ~£ = (xi, . . . ,x n ), 

2m 

and the corresponding (Hamiltonian) equations are x — §j3f and p — — f=r ■ In 
quantum mechanics, the momentum ~p* is replaced by = — zftV, the Hamiltonian 
operator by the Schrodinger (non-relativistic, stationary) operator 

and Hamilton's equations by the Schrodinger equation H*f> = . 
The eigcnfunction \P is known as the wave function, the eigenvalue E is the energy 
level, U(x) is the potential or potential energy. Solutions (^,E) of the Schrodinger 
equation are the states of the particle. In specific applications, the potentials will 
satisfy some conditions depending of physical restrictions such as barrier, scattering, 
etc., see [161 [27l [38l l42l [54] . 

The Schrodinger operator is known to be a self-adjoint operator: W = H when 
viewed within the context of a suitable complex and separable Hilbert space. As 
a consequence, the spectrum Spec(-ff) of H must be real, bounded below, and 
be the disjoint union of the point spectrum Spec p (ff) and the continuous spectrum 
Spec c (fl") (see for example [5l l50ll6T] ). Wave functions belonging to the continuous 
spectrum may arise from barrier potentials and from periodic boundary conditions. 
The transmission and reflection coefficients are related with the barrier potentials, 

In this paper, we consider the one-dimensional Schrodinger equation with the 
normalization h = 2m = 1, i.e. 

H^ = E^, H = -d 2 + V(z), 

where we may have z = x (cartesian coordinate) or z = r (radial coordinate). We 
let \t n denote the wave function corresponding to an energy level E = E n . 



Our purpose is to investigate solutions and their classification from the view- 
point of differential Galois theory. We should point out that this viewpoint has 
already been explored in nearby theories, for example by Spiridonov in [60] For 
sake of completeness, we now recall a few definitions related to quantum mechan- 
ics, and several useful results on the one-dimensional Schrodinger equation. 

Bound States A solution ^„ is called a bound state when its norm is finite 
and the corresponding eigenvalue E n belongs to the point spectrum of H , i.e. 

E n G Spec p (H), J \^ n (x)\ 2 dx < oo, n e Z + . 

Ground State and Excited States. Wave functions ^l, ■ ■ • , ^n, • • ■ of 

the bound states are numbered so that the corresponding energies are ordered 

Eq < Ei < ■ ■ ■ < E n < The state (^>o,Eq) with minimal energy is called the 

ground state and the remaining states ^i, . . . , . . . are called the excited states. 
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Scattering States. A solution ^ corresponding to the level energy E is called 
a scattering state when E belongs to the continuous spectrum of H and its norm 
is infinite. 

To model a particle moving in one dimension, we use the classical one dimen- 
sional Schrddinger equation with Cartesian coordinate x. We could also consider 
motion in three dimension by reducing the associated Schrddinger equation to the 
one-dimensional radial equation (see |27l 1381 154| ). 

Darboux Transformation. The Darboux transformation is a process which, 
starting from a second order linear differential equation such as the Schrddinger 
equation, produces a family of second order differential equations with a similar 
shape. This vague formulation is made clear in the following theorem, taken from 
from [20] , which is slightly more general than what is known today as the Darboux 
transformation. 

For any value of a constant m, let y = y m denote a solution of the equation 
dly + Pd x y + {Q- mR)y = 0. 
Suppose, in addition, that we are given a solution 8 of the equation 

dlO + Pd x 8 + Q0 = {). 

Then the function 



d x y-?fy 



is a solution of the equation 



Vr 

9VMJ \9VR 



dlu + Pd x u + ( 9VBB X [ — -= ) - eVRdl [ — != ) - mR ) u = 0, 



when m^O. 



This general Darboux transformation (y, 0) t— > u preserves the shape of the 
differential equation and maps Q to 9VRd x (jjj^j - 8y/Rd 2 x ■ In [20tl2T]. 

Darboux considered the particular case R = 1 and P = 0, which is known today as 
the Darboux transformation. It is the following corollary of the above: 

For any value of a constant m, let y = y m denote a solution of the equation 

dlv = (f( x ) + ™)V- 

Suppose, in addition, that we are given a solution 9 of the equation d^O = (f(x) + 
mi)0 for some fixed value m\ (e.g m\ — 0). Then the function 

u = d x y —y 

is a solution of the equation 

d 2 x u = Udl - mi + m) u, 
for m ^ mi . Furthermore, 

Ml (l)-mi = f{x) - 2d x ( ?f) = 2 ( MV - f(x) - 2m x . 
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Some examples can be found as exercises in the book of Ince |30l p. 132]. Darboux's 
transformation is related with with Delsarte's transformation operators (named 
transmutations in |22j ). which today are called intertwiners or intertwining oper- 
ators. A good short survey about Darboux transformations can be found in |52j . 



Crum, inspired by the works of Liouville 1391 140] obtained an iterative gener- 
alization of Darboux's result appropriate for Sturm-Liouvillc systems: he proved 
that the Sturm-Liouvillc conditions are preserved under Darboux transformations, 
see |19j . Crum's result is presented in the following theorem. Recall that the 
Wronskian determinant Wr of k functions f±, f%, . . . , /& is 

Wr(f u ...,f k ) = detA, A, s d\ : f r i,j = l,2,...,k. 

Crum's Theorem. Let ^2, • ■ • be solutions of the Schrodinger equation 
H 1 ^ = E~$> for energy levels E = E\ , E2, ■ ■ ■ , E n respectively. Letting 

9[n]= W £*^ '^*\ and V[n] = V - <ld% In Wr(* lt . . . , * n ), 

we obtain the Schrodinger equation 

H [n] V[n] = E^[n], where H [n] = -d\ + V[n] for E ^ E u 1 < i < n 

Crum's transformation coincides with Darboux's transformation when n = 1; 
as shown in [47] , iterations of Darboux transformations coincide with the Crum it- 
eration. Both formalisms allow us to obtain new families of Schrodinger equations 
while preserving the spectrum and the Sturm-Liouville conditions, see [421 147| . We 
note that there are extensions of Crum's iteration connecting the Sturm-Liouville 
theory with orthogonal polynomial theory |37j . 

Supersymmetric quantum mechanical system. Witten introduced, in 
[721 §6], some models for which dynamical breaking of super symmetry is possible. 

A supersymmetric quantum mechanical system is described by a Hamiltonian 
% together with several (self-adjoint) operators Qi which commute with the hamil- 
tonian ([QijH] = 0) and satisfy the anticommutation relations 

{Qi, Qj} = SijH, where {Qi, Qj} = QiQj + QjQi. 

The simplest example of a supersymmetric quantum mechanical system arises 
in the case n = 2, which is considered in this paper. The wave function of HQ = E$ 
is then a two-component Pauli spinor, 



The supercharges Qi are defined as 

o- 1 p±a 2 W(x) 
Q± = — 



= Qi, Q- =Q 2 , p = -id x , 



2 

where the superpotential W is an arbitrary function of x and o~i are the usual Pauli 
spin matrices. Using the expressions above we obtain %: 

H = 2Q^ = 2Ql = I ^ + I * W2 M + ^ W M , within (J J). 
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The supersymmetric partner Hamiltonians H± are given by 



l„o „ „ /M^\ 2 . 1 „ (W 



The potentials V± are called supersymmetric partner potentials and are linked 
with the superpotential W through a Riccati equation. Then T-L can be written as 



n = 



H+ 
H- 



which leads us to the Schrodingcr equations H + ^ + = E^ + and H_^>- = E^>^. 
To solve HQ = EQ is equivalent to solve simultaneously = E^+ and 

H-^- = 

In [721 §6], V.B. Matveev and M. Salle recast Darboux transformations in this 
context; they interpret the Darboux theorem as a Darboux covariance of a Sturm- 
Liouville problem and prove the following result (see also [521 §5-6]). 
The case n = 2 in Supersymetric Quantum Mechanics is equivalent to a single Dar- 
boux transformation. 

As a generalization of the method to solve the harmonic oscillator |28l 123] , the 
ladder (raising and lowering) operators are defined as 

A+ = -d x -^, A = d x -^, 

which are strongly related with the supercharges Q± in Witten's formalism. We 
have, 

A^ =0, A + A = AA + = H + = —d\ + V+(x), where 



V+(x) = V-(x) - 28 x = -V-(x) 



2 



*o ) V *o 

The supersymmetric partner potentials V+ and V— have the same energy levels, 
except for E$ ' = 0. In terms of the superpotential W{x), the operators A and A + 
are given by 

A + = -d x + W(x), A = d x + W(x). 

Similarly, the supersymmetric partner potentials V± (x) and the superpotential 
W(x) satisfy: 

V+(X)+ 2 V - {X) =W\X), [A,A+] = 2d x W(x). 

Shape invariance. Gcndcnshtc'in, in his remarkable paper |28j . introduced 
the concept of shape invariance, which is a property of families of potentials with 
respect to their parameters. Given a family a = (a n ) n of parameters, the shape 
invariance condition can be expressed as : 

V n +i(x; a n ) = V n (x;a„+i) + R(a n ), V- = V , V+ = Vi, 

where R is a remainder, which does not depends on x. When this relation holds, 
we say that V = V- = Vq is a shape invariant potential (with respect to the given 
family or transformation) , see also |52L I16|. 123] . 
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Gendenshtein Theorem. Consider the Schrddinger equation H^> = E n ^>, 
where V = V- = Vq is a shape invariant potential. If we fix the first level of energy 
Eq = ; then the excited spectrum and the wave functions are given respectively by 

n+1 n 

E n = ^i?(a fc ), *<-\x;ai) = Y[A + (x;a k )^ { ~ } (x;a n+1 ). 

k=2 k=l 

The structure of the rest of the paper is as follows. In section 2 we introduce the 
necessary concepts and results of the Differential Galois Theory of linear ordinary 
differential equations. In section 3 we introduce some terminology about solvability 
and we obtain some results about the invariance of the Galois group with respect 
to Darboux transformations. Section 4 is devoted to the study of the rational 
potentials and section 5 to transcendental potentials by means of the algebrization 
method. 

2. Differential Galois Theory 

In this section, we recall some differential Galois material to study the one- 
dimensional Schrbdinger equation, 

where V is some function and A is a constant (the spectral parameter). There 
are many references on differential Galois theory (or "Picard-Vessiot" theory), for 
example [TT I \18 [ 143 ] [65] . 

2.1. Picard-Vessiot Galois theory. We establish an algebraic model for 
functions and the corresponding Galois theory. 

Differential Fields. Let K be a commutative field of characteristic zero. 
A derivation of if is a map d x : K — > K satisfying d x (a + b) = d x a + d x b and 
d x (ab) = d x a -b + a- d x b for all a, b e K. We then say that (K, d x ) (or just K, when 
there is no ambiguity) is a differential field with the derivation d x . 
We assume that K contains an element x such that d x (x) = 1. Let C denote the 
field of constants of K: 

C = {ce K\d x c = 0}. 

It is also of characteristic zero and will be assumed to be algebraically closed. 
Throughout this paper, the coefficient field for a differential equation (e.g the 
Schrodinger equation d x ^f = (V — X)^) is defined as the smallest differential field 
containing all the coefficients of the equation. 

We will mostly analyze second order linear homogeneous differential equations, 
i.e equations of the form 

C := d x y + ad x y + by = 0, a, b £ K. 

so the rest of the theory will be explained in this context. 

Picard-Vessiot Extension. Let L be a differential field containing K (a 
differential extension of K). We say that L is a Picard-Vessiot extension of K 
for C if there exist two linearly independent yi , j/2 £ L solutions of £ such that 
L = K(yi,y 2 ) (i.e L = K(yi,y2,d x yi,d x y2)) and L and K has the same field of 
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constants C. 

In what follows, we choose a Picard-Vessiot extension and the term "solution 
of C will mean "solution of C in L" . So any solution of £ is a linear combination 
(over C) of 2/1 and y 2 - 

Differential Galois Groups A ^-automorphism a of the Picard-Vessiot ex- 
tension L is called a differential automorphism if it leaves K fixed and commutes 
with the derivation. This means that a(d x a) = d x (o~(a)) for all a € L and Va G K, 
a (a) = a. 

The group of all differential automorphisms of L over K is called the differential 
Galois group of L over K and is denoted by DG&l(L/K). 

Given a € DGal(£/-fT), we see that {ayi, cry 2 } are also solutions of C. Hence 
there exists a matrix 

4r=(° 5)eGL(2,C), 

such that 

2/2)) = (o-(yi) (7(2/2)) = (z/i 2/2) Ax- 

As a commutes with the derivation, this extends naturally to an action on a 
fundamental solution matrix of the companion first order system associated with 
C. 

,( yi v*\ = ( a (y^ \ = ( yi V2 \ A a 

XpxVi d x y 2 ) \o-(d x yi) o-(d x y 2 ) J [dxVi d x y 2 ) 

This defines a faithful representation DGal(L/if ) — > GL(2, C) and it is possible 
to consider DGal(L/i^) as a subgroup of GL(2, C). It depends on the choice of the 
fundamental system {2/1,2/2}, but only up to conjugacy. 

Recall that an algebraic group G is an algebraic manifold endowed with a group 
structure. Let GL(n,C) denote, as usual, the set of invertible n x n matrices with 
entries in C (and SL(n,C) be the set of matrices with determinant equal to 1). A 
linear algebraic group will be a subgroup of GL(n,C) equipped with a structure of 
algebraic group. One of the fundamental results of the Picard-Vessiot theory is the 
following theorem (see I33L 1351 ) . 

The differential Galois group DGal(L/i^) is an algebraic subgroup o/GL(2,C). 
In fact, the differential Galois group measures the algebraic relations between the 
solutions (and their derivatives). It is sometimes viewed as the object which should 
tell "what algebra sees of the dynamics of the solutions" . 

In an algebraic group G, the largest connected algebraic subgroup of G con- 
taining the identity, noted G° , is a normal subgroup of finite index. It is often 
called the connected component of the identity. If G = G° then G is a connected 
group. 

When G° satisfies some property, we say that G virtually satisfies this property. For 
example, virtually solvability of G means solvability of G° and virtual abelianity of 
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G means abelianity of G° (see |71j ). 

Lie-Kolchin Theorem. Let G C GL(2,C) be a virtually solvable group. Then 
G° is triangular izable, i.e it is conjugate to a subgroup of upper triangular matrices. 



2.2. Algebraic Subgroups of SL(2, C). We recall below some examples of 
subgroups of SL(2,C) for later use. 

Reducible subgroups These are the groups which leave a non-trivial sub- 
space of V invariant. They are classified in two categories. 

Diagonal groups: the identity group: {e} = | l) 1' ^ 6 n ~ roo ^ s ' = 

q c" 1 ) ' C " = ^ e ^-duplicative g rou P : G m = | c^ 1 ) ' C ^ ^ 
Triangular groups: the additive group; G a — i(n i)' ^ G w, the n — qu&si- 
roots: G{"> = jYg C" 1 ) ' C ™ = d e C }' thc Borel S roup: b = C ' kC = 



c 

c- 1 



Irreducible subgroups The infinite dihedral group (also called meta-abelian 
group): 

c G C* 1 U < [ I , d G C* i- and its finite subgroups 



v c' l J ' J \ °, 

D2„ (where c and c? spans the n-th roots of unity). 

There are also three other finite irreducible (primitive) groups: the tetrahcdral 
group A^ 2 of order 24, the octahedral group Sf L2 of order 48, and the icosahedral 
group Af La of order 120. 

Integrability. We say that the linear differential equation C is (Liouvillc) 
integrable if thc Picard-Vessiot extension L D K is obtained as a tower of differential 
fields K = Lq C L\ C • • • C L m = L such that L; = Lj_i(tj) for i = 1, . . . , m, 
where either 

(1) r] is algebraic over that is 77 satisfies a polynomial equation with 
coefficients in 

(2) 7] is primitive over L^-i, that is c^ry G J-h-\. 

(3) r/ is exponential over L,_x, that is d x r\jr\ £ Lj-i. 

We remark that the usual terminology in differential algebra for integrable 
equations is that the corresponding Picard-Vessiot extensions are called Liouvil- 
lian. The following theorem is due to Kolchin. 

The equation C is integrable if and only ifDGal(L/K) is virtually solvable. 



2.3. Eigenrings. Eigcnrings are a tool, introduced by M.F. Singer, to study 
the decomposability of systems or operators. A general reference is [6] (and refer- 
ences therein) . 
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Eigenrings of Systems. Consider the system d x X = —AX, denoted by [A], 
where A £ GL(2,K). The eigenring of the system [A], denoted by £(A), is the set 
of 2 x 2 matrices P in K satisfying 

d x P = PA- AP. 

Eigenrings of Operators. Let £ be a second order differential operator, i.e 
C := £(y) = 0. Denote V(£) as the solution space of C. The eigenring of £, 
denoted by £(£), is the set of all operators Dt of order less than that of £ for which 
JH(V(£)) is a subset of V(£), that is £3\ = 6£, where 6 is also an operator. For 
more details see [59l [6l [6HI [671 [68] . 

We study an equation £(y) = with £ = d x + pd x + q ; this is equivalent to 
the system [A] , where A is given by 

Relations between the system approach and the operator approach are given 
in the following simple lemmas . 

Lemma. Let £ = d 2 x + pd x + q, A = ( ° 1 ) and P = ( ° \) . Then: 



q p J \cd 

(1) If P e £(A), then m = a + bd x £ £(2). 

(2) If V\ = a + bd x ££(2), then P £ £(A), where P is given by 



P 



a b 

d x a — bq a + d x b — bp 



(3) 1 < dim c £(£) < 4 

~ b ^ f r- s 



(4) P£GL(2,K)^Zf -f +p^^-k q . 



Lemma. Assume that £ = d x + b, where b £ K. The following statements 
hold: 

(1) If dimc£(£) = 1, then either the differential Galois group is irreducible 
(Boo or primitive), or indecomposable (G CM, G non-diagonal). 

(2) If dime £(£) = 2, then either the differential Galois group is the additive 
group or it is contained in the multiplicative group (but is not the identity 
group). 

(3) If dime £(£) = 4, then the differential Galois group is the identity group. 
In this case we have 2 independent solutions £i and £2 in which £| 
and C1C2 are elements of the differential field K, i.e. the solutions of 
belong to to the differential field K. 

2.4. Some useful special functions. Some standard special functions will 
be used below. 

The Gauss hypergeometric equation. It is a particular case of Riemann's 
equation (see [32l SHI IHH SH [24]) and is given by 

2 (7 - (K + + l)g) a Kp 

d x y + 7z ; d x y -y = 0. 

x{\ — x) x{l — x) 

The cases when it is solvable can be found in the Kimura table I34J (or via Kovacic's 
algorithm, see appendix). 
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The confluent hypergeometric equations (Kummer, Whittaker). It 

is a degenerate form of the hypergeometric differential equation where two of the 
three regular singularities merge into an irregular singularity. For example, making 
"1 tend to oo" in a suitable way, the hypergeometric equation may transform to 
two classical forms: 

• Rummer's form 

_ 2 c — x _ a 

o x y + o x y y = 

x x 

• Whittaker 's form 

'1 k V 



d x y 



4 x Ax 2 



where the parameters of the two equations are linked by k = | — a and fx = | — -| . 
The Galoisian structure of these equations has been deeply studied in |41L 124) . The 
following theorem is due to Martinet & Ramis, [41] . 

The Whittaker 's differential equation is integrable if and only if either, k + /i€ 
\ + N. or k - n G \ + N, or — k + /i G \ + N, or —n - fx 6 \ + N. 

The Bessel equation. It is a particular case of the confluent hypergeometric 
equation and is given by 

„ 9 1 „ x 2 — n 2 

d 2 x y+-d x y+ — y = 0. 

x or 

Under a suitable transformation, the reduced form of the Bessel's equation is a 
particular case of the Whittaker's equation. Thus, we can obtain the following well 
known result, see [351 p. 417] and see also 



The Bessel's differential equation is integrable if and only if n g | + Z. 

By double confluence of the hypergeometric equation, that is making "0 and 1 
tend to oo" in a suitable way, one gets the parabolic cylinder equation (also known 
as Weber's equation): 

d x y = (7a; 2 - 7 - n ] y, 



which is integrable if and only if n G Z, see |36U24| . Under suitable transformations 
one can gets the Rehm's form of the Weber's equation: 

d x y = (ax 2 + 2bx + c)y, a^O, 
so that - — - is an odd integer. 

The hypergeometric equation, including confluences, is a particular case of the 
differential equation 

d 2 y+—d x y+—y, A G C, L = a Q + aix, Q = b a + b\x + b 2 x 2 . 

We mention this because the classical orthogonal polynomials and Bessel polynomi- 
als are solutions of such an equation (see |15 |, I31| . [46]). We remark that intcgrability 
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conditions and solutions of differential equations with solutions orthogonal polyno- 
mials, including Bcssel polynomials, can be obtained applying Kovacic's algorithm. 
Similarly, one can apply Kovacic's algorithm to obtain the same results given by 
Kimura |34j and Martinet & Ramis |41j . Also we recall that Duval & Loday- 
Richaud applied Kovacic's algorithm in |24j to some families of special functions, 
including several of the above. 

3. Notions of Algebraically Solvable and Quasi-Solvable Potentials 

The main object of our Galoisian analysis is the one-dimensional Schrodingcr 
equation, written as 

(3.1) C\ := = A*, H = -dl + V{x), V € K, 

where K is a differential field (with C as field of constants) . We are interested in the 
integrability (in Galoisian sense) of equation (|3.1j) . Denote by L\ the Picard-Vessiot 
extension of K for C\. The differential Galois group of C\ is DGal(L A /A). 

3.1. Algebraic Spectrum and Algebraic Solvability. Let A C C denote 
the set of eigenvalues A such that equation (|3.1[) is integrable in Galoisian sense. 
The set A will be called the algebraic spectrum (or alternatively the Liouvillian spec- 
tral set) of H. We remark that A can be 0, i.e., DGal(L A /A') = SL(2,C) VA G C. 
On the other hand, the classification of subgroups of SL2(C) shows that, if Ao G A, 
then (DGal(L Ao /A))° C B. 

We let A+ be the set {A G Anil : A > 0} and A_ be the set {A G AnM : A < 0}. 

Definition 3.1 (Algebraically Solvable and Quasi-Solvable Potentials). We say 
that the potential V(x) G K is: 

• an algebraically solvable potential when A is an infinite set, or 

• an algebraically guasi- solvable potential when A is a non-empty finite set, 
or 

• an algebraically non-solvable potential when A = 0. 

When Card(A) = 1, we say that V(x) G K is a trivial algebraically quasi-solvable 
potential. 

Examples. Let K = C(x). 

(1) If V(x) = x (Airy equation), then A = 0, V(x) is algebraically non- 
solvable, see [Ml[g6] , 

(2) If V(x) = 0, then A = C, i.e., V(x) is algebraically solvable. Furthermore, 

DGal(L /A) = {e}, and DGal(L A /A) = G m for A ^ 0. 

(3) If V(x) = ^- + |, then A = {n : n G Z}, V(x) is algebraically solvable. 
This example corresponds to the Weber's equation given previously. 

(4) If V(x) = x 4 — 2x, then A = {0} and V{x) is a trivial algebraically quasi- 
solvable potential (see next section for tools to prove this). 

We are interested in the spectrum (analytic spectrum) of the algebraically solvable 
and quasi-solvable potentials, that is, Spec(iJ) flA^B. 
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Potentials for which Spec(if)(~lA is an infinite set are called solvable (or exactly 
solvable) potentials (see Natanzon [45] ) in the physics litterature. Similarly, poten- 
tials for which Spec (H) n A is a finite set are usually called quasi- exactly solvable (or 
quasi-solvable) potentials (Turbiner [63] . Bender & Dunne [§], Bender & Boettchcr 
[7], Saad et al. [53], Gibbons & Vesselov [29]). 

3.2. Changing the Base Field. In what follows, we will make transforma- 
tions which may change the base field; typically, when K is not C(x), we will look 
for transformations to another differential with rational coefficients. We now collect 
facts on how this affects the differential Galois group and its solvability. 

Definition 3.2. Let £ and £ be a pair of linear differential equations defined over 
differential fields K and K respectively, with Picard-Vessiot extensions L and L. 
Let tp be a transformation mapping K to K and £ to £ (so we may also choose tp 
to map L to IS). We say that: 

(1) tp is an isogaloisian transformation if 

DGal(L/A") = DGal(L/K). 

When L = I and K = K, we'll say that tp is a strong isogaloisian trans- 
formation. 

(2) tp is a virtually isogaloisian transformation if 

(DG&\(L/K))° = (DGal(Z/#))°. 

Remark 3.3. The eigenring s of two operators £ and £ are preserved under iso- 
galoisian transformations. 

Proposition 3.4. Consider the differential equations 

£ := d 2 x y + ad x y + by = 0, £ := d*( = r(, a, b, r G K. 

let k G Q, / G -ft', a = 2ft<9 z .(ln/) and tp be the transformation such that C h- > C. 
The following statements holds: 

(1) tp is a strong isogaloisian transformation for c£Z. 

(2) tp is a virtually strong isogaloisian transformation for kg Q \ Z. 

Proof. Let B = {2/1,2/2} be a basis of solutions of C and L be the corre- 
sponding Picard-Vessiot extension of K; similarly, let B' = { Ci 3 C2 } be a basis of 
solutions of £ if a Picard-Vessiot extension L. With the change of dependent vari- 
able 2/ = C e ~ 3 ^°we obtain r = a 2 /4+d x a/2 — b and K = K. Thus, the relationship 
between L and L depends on a: 

(1) If k = n G Z, then = {.f n yi, f n y2} which means that I = L and 93 is 
strong isogaloisian. 

(2) If « = with gcd(n,m) = 1, £ £ Z, then 6' = {f^yi,f^y 2 } which 
means that £ is either an algebraic extension of degree at most m of I, 
and y> is virtually strong isogaloisian, or I = I when f™ £ K which 
means that tp is strong isogaloisian. 

□ 



GALOISIAN APPROACH TO INTEGRABILITY OF SCHRODINGER EQUATION 13 



Remark 3.5. The transformation (p in proposition 13.41 is not injective, there are a 
lot of differential equations C that are transformed in the same differential equation 
£. They are called projectively equivalent. 

The following is an easy consequence of the proposition: 

Corollary 3.6. Let C be the differential equation 

d x (ad x y) = (Xb - ff)y, a, & £ K, A,/ieC 



where L, L, C and <p are given as in proposition \3.4\ Then either L is a quadratic 
extension of L ( which means that if is virtually strong isogaloisian ) or L = L ( this 
is when a' £ K) which means that Lp is strong isogaloisian. 

4. Galoisian Aspects of Darboux Transformations and Shape Invariance 

Here we present a Galoisian approach to the Darboux transformation, the Crum 
iteration and shape invariant potentials. 

Recall that Wr(yi, . . . , y n ) denotes the Wronskian 



Wr(y 1 , ...,y n ) 



yi ■■■ yn 



£)n—l n . Qn-L, 
"x yi ■ d x Vn 

The Darboux transformation will be abreviated by DT and its n-th iteration by 
DT n ; similarly, the abreviation CI n will refer to the Crum iteration. 

As above, K denotes the coefficient field of the studied linear differential equa- 
tion, i.e. the smallest differential containing the coefficients of the equation. We 
mostly consider equations over the base field K = C(x). 

4.1. Darboux Transformations. The following result comes from the origi- 
nal Darboux transformation (see 1201 121| ). but rewritten in a way which is adapted 
to our Galoisian analysis. 

Theorem 4.1 (Darboux Transformation). Consider the Schrddinger operators 
H± = &x + y±{ x )i with V±(x) £ K, and assume that the spectrum satisfies 
A 0. Let C\ (resp. C\) denote the Schrddinger equation H-^f^ = A\I/( _ ) 
(resp. i2+*(+) = A*(+) ). 

Let DT be the transformation such that C H> C, V- M> V+, , 3 /( ~ ) M> \I>( + ). Then the 
following statements holds: 

i) : DT(F_) = V + = <^ (i) +\ 1 = V-- 2^0* *<">), 

DT^^ ') = = —^=j, where Vf^ ' is a particular solution of C\ 1 , 

1 1 *Ai 1 

Ai e A. 

ii) : DT(*(")) = *(+) = x *<-> - ^(ln = ^g^ 2 , A ? \ l} 

where v]/( _ ) = ty{ ' is the general solution of £\ for A £ A \ {Ai} and 
^r(+) = tyy' is the general solution of C\ also for A £ A \ {Ai}. 

PROOF. We treat each item separately. 
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i): By hypothesis, we have Card(A) > 1, Ai € A, 

.ff_*^ = Ai*^~\ and setting H± = H± — Ai, V± = V± - Ai we obtain 

H-^ x ^ — 0. Thus, the superpotential W = — d x ]n^ x , depending of 
Ai, exists and satisfies 

(-a x + w) (d x + w) * j-> = = o. 

Now, developing both sides of the expression, we have V- = — d x W + W 2 . 
So, interchanging operators we have 

(a, + w) (-a x + w) = = o, 

so that developing both sides of the expression we have V+ = d x W + W 2 

'it 



which lead us to W — d x ln\l/( + ' ) . Thus, we obtain 



*(-)" 

Ai 



On another hand, as V+ - V_ = V+ - V_ = 2d x W, V + + V_ = V + + V_ - 
2Ai = 2W 2 , we obtain DT(V_) = V + which is given by 

ii): Suppose that Card(A) > 1, A ^ Ai and *( ± ) = ^[ ±} satisfying 

= A* (_) , A^Ai, 
applying again the raising and lowering operators we have 

{-d x + W) (d x + W) <~> ? 0, (d x + W) {-d x + W) * (+) ^ 0, 
so that 

(-d x + w) {d x + w) *(-) = {-d x + w) ¥+\ 

we obtain 

*(+) = [d x + W) * (_) = 8 x tf (_) + 
= a vjj(-) Al ^r(-) = v Al ' L 

Now, we will show that effectively * (+) satisfy = A* (+) for 

A ^ Ai. Let suppose that 

*(+)=^(-)-*(-)^^, 
< } 

we can see that 

fl£*(-> _ = $(-)( Al _ A ). 
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therefore, differentiating^ ( + \ we have 



By differentiating d x ^ + >, we have 



(-) \ 2 Q ,T f (-) «2,T,(-) 



d 2 x ¥ +) = d x ^-\X! - A) + 20 x ¥^ ^ +2*(") Al x x 



\ Al / Ai Al 



\ Ai / \ Ai / Ai 

a x *(-)(Ai-A)+29^-) I ^^-j -2*(-)|^ T ^-j _fl x tf(-)JL_^_ 

V Ai / V ^Ai / Ai 

x *C->_*<-)^* ] 




+ Ai — A d x ¥^ - 



~ -^TT" + Ai - A | 

y_ + 2Ai - aJ * (+) , 

which implies that, for A ^ Ai, the function v[r(+) satisfies the Schrodinger 
equation = A* (+) , 
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□ 



Remark 4.2. According to the corollary of the general Darboux transformation, 
see section[TJ we have m = —A, mi = — Ai, f{x) = V-,y = and u = y = . 
This gives us a faster proof of theorem 14.11 

In agreement with the previous theorem, we obtain the following results. 

Proposition 4.3. DT is isogaloisian and virtually strong isogaloisian. Further- 
more, if d x (hi^ Xi ) £ K, then DT is strong isogaloisian. 

Proof. Let K and A A be respectively the coefficient field and the Picard- 
Vessiot extension of the equation C\. Similarly, let A, L\ be the coefficient field 
and the Picard-Vessiot extension of the equation C\. As DT(V_) = V + = 2W 2 — 

V- - 2Ai, where W = -^(ln*^), we have K = K (d x (ln V^)} ■ The Riccati 
equation d x W = V- — W 2 has one algebraic solution (Singer 1981, |57j ). in this 
case W = —dxQja.^^). Let be a basis of solutions for equation C\ 

and (\f , ;^», ^0^0 a basis of solutions for equation C\. Since the coefficient field 
for equation C\ is K = K (d x (ln^^)) , we have L = K{#t\)>^(2 \)) and 

i=^si ) ^ai))=Mti)'*si)'^( in *i7 ) )) 

for A = Ai and for A ^ Ai. Since ^(lnVP^) is algebraic over A, then 

(DGal(A A /A))° = (DGal(Z A /A))°, DGal(A A /A) = DGal(Z A /A), 

which means that DT is a virtually strong and isogaloisian transformation. 

In the case d x (ln^^') € K, then A' = K and A = L, which means that DT is a 

strong isogaloisian transformation. □ 

Proposition 4.4. Consider Schrodinger operators £ A := H- — A and £ A := H+ — A 
such that DT(iA_ — A) = A + — A. The eigenrings o/£ A and £ A are isomorphic. 

PROOF. Let £(£\) and £(£\) denote the eigenrings of £ A and £ A respectively. 
By proposition 14.31 the connected identity component of the Galois group is pre- 
served by Darboux transformation. Let T <G £(£\). Let Sol(£ A ) and Sol(£ A ) be the 
solution spaces for £\ty(~*> = and £^ = respectively. To transform £(£\) 
into £(£ A ), we follow the commuting diagram: 



Sol(£ A ) 

At 

Sol(£ A ) 



Sol(£ A ) 

A 

Sol(£ A ) 



so T = ATA^ mod(£ A ) e £ (£ A ) , 



where A^ and A are the raising and lowering operators. 



□ 
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Example. Consider the Schrodinger equation C\ with potential V = V- = 0, 
which means that A = C. If we choose Ai = and as particular solution \1/q = x, 
then for A 7^ the general solution is given by 

Applying the Darboux transformation DT, we have DT(£ A ) = £\, where 
and for A 7^ 



DTf*") = = ci(V^Xx-l)e^ _ c 2 (y^Ax+l) e -^ 

A A X X 

We can see that K — K = C(x) for all A £ A and the Picard-Vessiot exten- 
sions are given by L = L Q = C{x), L x = L\ = C(x,e^ x ) for A £ C*. It 
foUowsthat DGal(L /K) = DGal(Z /A') = e; for A ^ 0, we have DGal(L A /A') = 
r>Gal(L\/ K) = G m . The eigenrings of the operators £0 an d £a are given by 

£(£0) = Vect (1, d x , xd x - 1, x 2 d x - x) , 



£(£ ) =Vcct ( l,xd x -l,x 4 d x -2x 3 ,^ + A 



X' X" 



and for A 7^ 



£(£ A ) = Vect (1, d x ) , £(£ X ) = Vect (l, - ^A + A J ^ _ A 
where £ A := £a* ( ~ ) = and C\ := £x* (+) = 0. 

Applying iteratively the Darboux transformation, theorem 14. 1 L and by propo- 
sitions 231 S3] we have the following results. 

Proposition 4.5 (Galoisian version of DT„). We continue with the notations of 
theorem\4J\ We set V = VL, i? (0) = H-, * (0) = *(~) We still assume that 
the algebraic spectrum A is non-empty. We consider the operators (given 
by F( n )*(") = AtfWj obtained by successive Darboux transformations. We let 
K n be the coefficient field of (i.e V n € K n ) and K = Kq . Let denote 
the Picard- Vessiot extension of . The Darboux transformation DT„ is such 
that £ A n) i-> £ A ™ +1 \ V n o Vn+i, * A n) h> * A n+1) and L A " +1) the Picard- Vessiot 

extension of £ A . 

Then the following statements hold: 

i) : DT„(F_) = DT(F„) = K+i = (ln*i n) ) = K_-2£Lo^ (l**!?), 

(AO 

where w A is a particular solution for A = Afe, A; = 0, . . . , n. /n particular, 

if A„ = Ao afid A = C, £/ien f/iere exists ^\ such that V n 7^ V r „„2 ; wii/i 
n > 2. 

ii) : DT(*W) = DT„(^-)) = ^ = flt*W-*W^ = 
where is a general solution for A € A \ {A„} of 
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Hi): K n+1 =K n (d x {bx^)y 

iv) : DT„ is isogaloisian and virtually strongly isogaloisian. Furthermore, if 
9 x (ln\&^ ) € K n then DT„ is strongly isogaloisian. 

v) : The eigenrings of — A and — A are isomorphic. 



Proof. By induction on theorem 14.11 wc obtain i) and ii). By induction on 
proposition 14.31 we obtain iii) and iv). By induction on proposition 14.41 we obtain 
v). 

□ 



Example. Starting with V = 0, the following potentials can be obtained using 
Darboux iteration DT„ (see [iUl,[i~2]). 

/) v n = n ( n ~ ^ //) Vn - m2n ( n ~ !)( &2 ~ « 2 ) 



(bx + c) 2 ' " (acosh(?7ia;) + 6sinh(ma;)) 2 ' 

Aabm 2 n(n — 1) . m 2 n(n — l)(b 2 + a 2 ) 



IH)Vn= ,__„,, , IV) V, 



(ae mx + be mx ) 2 ' (ocos(mi) + 6sin(ma;)) 2 

In particular for the rational potential given in /), we have K = K n = <C(x) and 
for A„ = A = 0, we have 

*o l) = Cl + c 2(bx + c) n+1 , so that DGal(L /A") = DGal(4 n) /A) = e, 
(ox I c) 

£{HW) = Vcct (l, xd x 1, x 2n+2 d x (n + l)x 2n +\ ^ + -^j , 

whilst for A ^ and A„ = 0, the general solution is given by 

*X ( x ) = cif n (x, \)h n (sm(<f\x) + c 2 g„{x, X)j n (cos(\/Xx), 
where f n ,g n ,h n ,j n £ C(x), so that 

DG&1(L\/K) = DG&l(L^ l) /K) = G m , 

and 

dime S(H - A) = dime E{H {n) - A) = 2. 

4.2. Crum Iteration. 

Proposition 4.6 (Galoisian version of CI„). Consider C\ given by = A^, 
H = -d 2 + V, V G K, such that Card(A) > n for a fixed n e Z+. Lei £ A ™ } be 
given by = A*( n ), w/iere ffW = d 2 x + V„, V n € K n . Let CI„ &e fee 

transformation such that C\ h-> £ A , V h-> Ki, ($An • ■ • , ^A n) *a) ^ ^ A > where 
for k = 1, . . . , n and i/ie equation C\, the function ty\ is the general solution for 
A 7^ Afc and \& Aj! *s a particular solution for A = A^. TTien i/ie following statements 
holds: 

i) : CI„(£ A ) = £i n) w/iere CI„(y) = V„ = V-2d 2 . . . , * A J) and 

CIn( ^ A) ~* A - W(* Al> ...,*,J ' 

where * s ^ e general solution of ■ 

ii) : K n = K(d x (In W(tt Al , . . . , tt A J)) . 
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iii) : CI„ is isogaloisian and virtually strongly isogaloisian. Furthermore, if 

ft,(lnW(¥A 1 ,... > tt An )) 6 Kn, 

then Gl n is strongly isogaloisian. 

iv) : The eigenrings of H — A and H^ n > — A are isomorphic. 

PROOF. By induction on theorem l4.1l we obtain i). By induction on proposition 
we obtain ii) and iii). By induction on proposition 14.41 we obtain iv). 

□ 



Example. To illustrate the Crum iteration with rational potentials, we consider 
V = Jr. The general solution of C\ := H^> = A* is 

Cxe kx {kx - 1) c 2 e~ kx (kx + 1) 2 

T" ; A = K , 

X X 

the eigenfunctions for Ai = — 1, and A2 = —4, are respectively given by 

e- x (x + l) _ e- 2x (2x + l) 

x 2x 

Thus, we obtain 

Ch(V) = V 2 = 



(2x + 3) 2 

and the general solution of £ A 2) := i?( 2 )*( 2 ) = A* (2) is 
PT ftt ^ - *( 2 ) - Cl (M2* + 3)-2)e fc * c 2 (2 + fc(2x + 3))e-^ _ 2 

The differential Galois groups and eigenrings are given by: 

DGal(L //\) = DGal(-4 2) /-ftT) = e, dim c £(if) = dim c £(# (2) ) = 4, 
and for A 7^ 

DGel(L\/K) = DGal(L^ / K) = G m , dim c £(H - A) = dim c £(ii^ - A) = 2. 

Proposition 4.7. TTie super symmetric partner potentials V± are rational functions 
if and only if the superpotential W is a rational function. 

PROOF. The supersymmetric partner potentials V± are written as V± = W 2 ± 
d x W. We start considering the superpotential W G C(x), so trivially we have 
V± G C(x). Now assuming that V± G C(x) we have d x W G C(x) and Vl^ 2 G C(x), 
which implies that Wd x W G C(x) and therefore W G C(x). 

□ 



Corollary 4.8. The superpotential W G C(x) if and only if DT is strong isoga- 
loisian. 

Proof. Assume that the superpotential W G C(x). Thus, by proposition 14.31 
DT is strong isogaloisian. Now, assume that DT is strong isogaloisian. Thus, 
V± G C(x) and by proposition 14.71 we have W G C(x). 

□ 
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4.3. Shape Invariance. The following definition is a partial Galoisian adap- 
tation of the original definition given in |28j (K = C(x)). The complete Galoisian 
adaptation is given when K is any differential field. 

Definition 4.9 (Rational Shape invariant Potentials). Assume V±(x; n) £ C(x; //), 
where fx is a family of parameters. The potential V = V- £ C(x) is called a rational 
shape invariant potential with respect to \x and E = E n if there exists / such that 

n+l 

V+(x;a ) = V-(x;ai) + R(ai), a 1 = f(a ), E n ^^R(a k ), E =0. 

k=2 

Remark 4.10. We propose the following steps to check whether V £ C(x) is shape 
invariant. 

Step 1.: Introduce parameters in W(x) to obtain W(x; fj,), write V± (x; fj,) = 

W 2 (x; n) ± d x W(x; //), and replace \i by ao and oi. 
Step 2.: Obtain polynomials V £ C[x; ao, ai] and Q £ C[x; ao, ai] such that 

d x (V + (x;a ) - V-(x;ax)) = — r. 

Q(a;;ao,ai) 

Step 3.: Set V(x; ao, ai) = 0, as polynomial in x, to obtain a\ in function 
of ao, i.e., ai = /(ao). Also obtain R(a{) = V+(x;ao) — V-(x;a\) and 
verify that exists k e Z+ such that R(ai) + • • • + R(dk) ^ 0. 

Example. Consider the superpotential of the three dimensional harmonic oscillator 
W(r;£) = r — ^tl. By step 1, the supersymmetric partner potentials are 

V-(r;l)=r> + e -^l-2l-3, V + (r;i) = r> + « + J + ^ - 21- 1. 

By step 2, we have d r (V+{r;a ) - V-(r;ai)) = -2 a ° +3a %a'~ Ql+2 - By step 3, 
(a + l)(a + 2) = 01(01 + 1), so that ai = f(a ) =a + l, a n = /(o n _i) = a + n, 
i?(ai) = 2. Thus, we obtain the energy levels E n = 2n and the wave functions 
$(~\r;£) = (r;£) ■ ■ ■ A^(r;£ + n - 1)* _) (r;t + n), compare with [23] , 

By theorem 14. H and propositions 14. 31 14.41 and 14.71 we have the following result. 

Theorem 4.11. Consider C n := = E n ^^' with Picard-Vessiot extension 

L n , where n G Z+. If V = V— € C(x) is a shape invariant potential with respect to 
E = E n , then 

DGal(L n+1 /K) = DGal(L n /K), £{H - E n+1 ) ~ £{H - E n ), n > 0. 

Remark 4.12. As W £ C(x), a differential automorphism a commutes with the 
raising and lowering operators A and A^ . Furthermore the wave functions ^1 ' 
can be written as ' = P n f n il! , where P n is a polynomial of degree n in x 
and f n is a sequence of functions with fo(x) = 1 such as was shown in the case of 
Harmonic oscillators and Coulomb potentials. 
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5. Schrodinger Equations with Rational Potentials 

Along this section we take K = C(x) to be our coefficient field. So all our 
differential operators will have coefficients in K = C(x). 

5.1. Polynomial Potentials. We start by considering the Schrodinger equa- 
tion (|3.1[) with polynomial potentials, i.e. V £ C[x], see |13l 1641 . For simplicity 
and without loss of generality, we consider monic polynomials (if the polynomial is 
c„x™ + . . . + c\x + Co, apply the transformation x i-> n+ ^fc~^x). 

When V is a polynomial of odd degree, is well known that the differential Ga- 
lois group of such a Schrodinger equation (|3.1[) is SL(2, C), see 



We present here the complete result for the Schrodinger equation (|3.1[) with 
non-constant polynomial potential (Theorem 15. 2[) . see also [H §2]. The following 
lemma is useful for our purposes. 

Lemma 5.1 (Completing Squares, [4]). Every even degree monic polynomial of 
can be written in one only way completing squares, that is, 

2n-l / n-1 \ 2 n-1 

(5.1) Q 2 „(x) = x 2n + <lk* k = + £ ^ + bkxk > 

k=0 \ k=0 ) k=0 

where 

Q.2n-l _ (fen-2 — a n-l _ Q2n-3 — 2a„_ia„_ 2 

O-n-1 — 2 ' a "-2 — 5 ' a ™-3 — 2 ' ' ' ' ' 

q n — 2aia n _i — 2a2a„_2 — • • • , 2 l n 

a = , = 9o-ao' 0i = gi-2a ai, 

Ori-i = 9n-i — 2aoa„„i — 2aia„_2 — • • • . 



PROOF. See lemma 2.4 in [H p. 275]. 

□ 

We remark that V(x) as in equation (|5.ip can be written in terms of the su- 
perpotential W(x), i.e., V(x) = W 2 (x) — d x W(x), when 

n— 1 n— 1 

fc=l fc=0 

and W(x) is given by 

n-l 

x" + a fc x fc . 

fe=0 

The following theorem also can be found in [4) §2], see also [3]. Here we present 
a quantum mechanics adapted version. 



Theorem 5.2 (Polynomial potentials and Galois groups, [4]). Let us consider the 
Schrodinger equation (|3.ip , mf/i V(x) € C[x] a polynomial of degree k > 0. Then, 
its differential Galois group DGal(L\/ ' K) falls in one of the following cases: 
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(1) DGal(L x /K) = SL(2,C), 

(2) DGal(L x /K) = B, 

and the eigenring of H — X is trivial, i.e., £(H — A) = Vect(l). Furthermore, 
DGal(L\/ K) = B if and only if the following conditions hold: 

(1) V(x) — X is a polynomial of degree k = 2n writing in the form of equation 

USD. 

(2) &„_i — n or — b n —i — n is a positive even number 2m, m G Z+. 

(3) There exist a monic polynomial P m of degree m, satisfying 

/ n-l \ / n-2 n-l \ 

dlPm + 2 I x n + a fc x fe diP™ + na-™ _1 + ^(fc + l)a fc+ ix fc - ^ fe fe x fe P m = 0, 

\ fe=0 / V k=0 k=0 I 

or 

I n-l \ / n-2 n-l \ 

9^P m - 2 i" + ^ a k x k d x P m - nx"^ 1 + y^(fc + l)a fc+ ix fe + ^ & fc a; fc P m = 0. 

V fe=0 / V k=0 k=0 ) 

In such cases, the only possibilities for eigenfunctions with rational superpotentials 
are given by 



* A =P m e^, or * x = P m e- f <>*\ where f(x) = — — + V 

n + 1 f—i 



n— 1 in II 



fe =o fc + 1 



PROOF. See theorem 2.5 in gj p. 276]. 

□ 

An easy consequence of the above theorem is the following. 



Corollary 5.3. Assume that V(x) is an algebraically solvable polynomial potential. 
Then V(x) is of degree 2. 

Proof. Writing V(x) — X in the form of equation (15.1[) wc see that — n = 
2m or — b n —i — n = 2m, where m G Z + . Thus, the integrability of the Schrodinger 
equation with Card(A) > 1 is obtained when 6„_i is constant, so n = 1. □ 



Remark 5.4. Given a polynomial potential V(x) such that Spec p (iJ) H A =^ 0, wc 
can obtain bound states and normalized wave functions if and only if the potential 
V(x) is a polynomial of degree 4n + 2. Furthermore, one integrability condition 
of H^f = A 'J for A e A is that b% n must be an odd integer. In particular, if the 
potential 

V{x) =x 4n+2n +fix 2n , n >0 

is a quasi-exactly solvable, then [i is an odd integer. For this kind of potentials, we 
obtain bound states only when [i is a negative odd integer. 

On another hand, the non-constant polynomial potentials V(x) of degree 4n 
are associated to non-hermitian Hamiltonians and VT invariance which are not 
considered here, see [JJ. Furthermore, one integrability condition of = X^> for 
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A G A is that &2n-i must be an even integer. In particular, if the Schrodinger 
equation 

H^ = X^, V(x) = x 4n + f ix 2n -\ A G A 
is integrable, then \i is an even integer. 

We present the following examples to illustrate the previous theorem and re- 
mark. 



Weber's Equation and Harmonic Oscillator. The Schrodinger equation with 
potential V(x) = x 2 + q\x + qa corresponds to the Rehm form of the Weber's 
equation. By lemma [57T1 we have 

V(x) - A = {x + a ) 2 + b , ao = qx/2, b = q - q\/A- X. 

So that we obtain ±60 — 1 = 2m, where m G Z+. If 60 is an odd integer, then 

~DGal(L x /K) = B, £(H - X) = Vcct(l), A G A = {±(2m+l) + q -q 2 /4 : m G Z + } 

and the set of eigenfunctions is either 

*a = P m e^ x2+qiX \ or, * A = P m e~^ + ^ x \ 

In the second case we have bound state wave function and Spec p (i?)nA = Spec p (iJ) = 
{E m — 2m + 1 + go — Si/4 : m G Z+}, which is infinite. The polynomials P m are 
related with the Hermite polynomials H m , I15L I31L 1461 . 

In particular we have the harmonic oscillator potential 

V(x) = \u 2 x 2 |, 

and the Schrodinger equation H*S> = E^f. Through the change of independent 

variable x 1— > \J~£;x we obtain V(x) = x 2 — 1 and A = ~E, that is, q\ = and 

50 = — 1. It follows that A = {±(2m+l) — 1 : m G Z + } and the set of eigenfunctions 
is either 

*a = P m e^ x \ or, * A = P m e~^ x \ 
where as below, DGal(L A /A') = B and £{H - A) = {1} for all A G A. In the second 
case we have bound state wave function, Spec p (7J) n A = Spec p (i?) = A + = {2m : 
m G Z + } and P m = H m . The wave functions of H^> = E^> for the harmonic 
oscillator potential are given by 

*m = H m I \l-x) $0, *o = e"^ 2 , E = E m = muj. 




Quartic and Sextic Anharmonic Oscillator. The Schrodinger equation with 
potential V(x) = x i +q3X 3 +q2X 2 +q±x+qo can be obtained through transformations 
of confluent Heun's equation, though we will not adress this here. By lemma 15.11 
we haw 

V(x) — A = (x 2 + a\x + ciq) 2 + b\x + bo, 

where a\ = q^/2, ao = (72/2 — a\/2, b\ = q\ — 2a^a\ and 60 = <7o — Oq — A. So we 
obtain ±b\ — 2 = 2m, where m G Z + . If A ^ 0, then 61 is an even integer, P m satisfy 
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the relation (jA.lj) and DGal(L\/ K) = B for all A e A. The set of eigenfunctions is 
either 

* A - P m e-+— +a " x , or, * A = P m e \ J, 

where A and m are related, which means that A is finite, i.e., the potential is alge- 
braically quasi-solvable. In particular for 33 = 2il, q2 = I 2 — 2k, q\ = 2i(lk — J) and 
go = 0, we have the quartic anharmonic oscillator potential, which can be found in 

m- 

Now, considering the potentials V(x, fi) = x 4 +4x 3 + 2x 2 — fix, again by lemma 
15.11 we have 

V(x, fi) - A = {x 2 + 2x - l) 2 + (4 - fi)x - 1 - A, 
so that ±(4 — fi) — 2 = 2n, where n € Z + and in consequence fi £ 21. Such fi can 
be either fi = 2 — 2n or fi = 2n + 6, where n £ Z + . By theorem 15.21 there exists a 
monic polynomial P n satisfying respectively 

d 2 x P n + (2x 2 +4x- 2)d x P n + ((fi - 2)x + 3 + X)P n =0, fi = 2 - 2n, or 

d 2 x P n - (2x 2 + 4x - 2)d x P n + ((fi - 6)x - 1 + A)P n =0, fi = 2n + 6 
for A 7^ 0. This algebraic relation between the coefficients of P n , fi and A give us 
the set A in the following way: 

(1) Write P n = x n + c n -\x n ~ l + . . . + cq, where Cj are unknown. 

(2) Pick fi and replace P n in the algebraic relation (|A.1[) to obtain a polyno- 
mial of degree n with n+1 undetermined coefficients involving c , . . . c„_i 
and A. Each of such coefficients must be zero. 

(3) The term n + 1 is linear in A and c n _i, thus we write Cn—i in terms of A. 
After of the elimination of the term n + 1, we replace c„_i in the term n to 
obtain a quadratic polynomial in A and so on until we reach the constant 
term which is a polynomial Q„ + i(A) of degree n + 1 in A. It follows that 
A = {A : Q„+i(A) = 0} and Co, . . . , c„_i are determined for each value of 
A. 

For fi = 2n + 6, we have: 

n = 0, V(x,6), P Q = 1, A = {1} 

n=l, V(a;,8), Pi=x + 1tV2, A = {3 ±2^2} 

and the set of eigenfunctions is 

*A,p = P n e-^ 3 ~ x2+x . 

In the same way, we can obtain A, P n and ^a.^ for fi = 2 — 2n. However, we do 
not have bound states, Spec p (#)nA = 0, DG&\(L\/K) = B and £(H—\) = Vect(l) 
for all A G A. 

The well known sextic anharmonic oscillator x e + q$x 5 + • • • + q%x + cjq can be 
treated in a similar way, obtaining bound states wave functions and the Bender- 
Dunne orthogonal polynomials, which corresponds to Q n +i(A), i.e. we also recover 
the results of [HI 1291 153] . The Schrodingcr equation with this potential, under 
suitable transformations, also falls in the class of confluent Hcun equations. 
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5.2. Rational Potentials and Kovacic's Algorithm. Here we apply Ko- 
vacic's algorithm (see appendix A) to solve the Schrodinger equation with rational 
(non-polynomial) potentials. 

Three dimensional harmonic oscillator potential: 

!. 2 2 , + ( t , 3 



In this case, the Schrodinger equation is written as 



By the change r4 (^y^Jrwe obtain the Schrodinger equation 

%* = (r2 + W±*l-(2l + 3)-\\% X=-E 
\ r 2 J to 

and, in order to apply Kovacic's algorithm, we set 

R = r 2 + £ Jl±ll _ (2^ + 3)- A. 



We can see that this equation could fall in case 1, in case 2 or in case 4 of 
Kovacic's algorithm. We start discarding the case 2 because by step 1 (of Kovacic's 
algorithm) we should have conditions C2 and 003; so wc should have E c = {2,4 + 
4:£,—4:£} and Eoo — {—2}, and by step 2, wc would have n = —4 ^ Z+, so that 
D = 0. So this Schrodinger equation never falls in case 2. Now, we only work with 
case 1; by step 1, conditions C2 and 003 are satisfied, so that 

± l±(2l+l) r/-i ± T (A + 2l + 3)-l 
= °. a ? = o ' =r, a* = . 

By step 2 we have the following possibilities for n £ Z + and for A G A: 



A++) 


n 


= a to 


-<*o = 




+ 6 + A), 


A = -In -4^-6 


A+-) 


n 


= a to 




-|(4H 


-A), 


A = -2n-4, 


A-+) 


n 


= a x 


-at = 


A 
2 ' 




A = 2n, 


A-) 


n 


= a^o 




i(4£ + 


2 + A), 


A = 2n - 4£ - 2, 



where A ++ U A H U A v U A = A, which means that A = 2m, vn 6 Z. Now, for 

A 6 A. the rational function oj in Kovacic's algorithm is given by: 

A++) c = r + ^i, 

A+_) w = r-$, 

A_+) c = -r+^±, 

A—) o; = -r-£, 





= r 2 A 


t(Z+l) 


+ (21 + 


3) 


+ 2n, 




— r 2 _| 


1(1+1) 


- (21- 


1) 


+ 2n, 


Rn 


= r 2 H 


1(1+1) 


- (2£ + 


3) 


- 2n, 


Rn 


= r 2 H 


1(1+1) 


+ (2l- 


1) 


- 2n, 
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where R n is the coefficient of the differential equation C n := d 2 ^> = R n ^, which 
is integrable for every n. For every A G A, we can see that DGal(L„/ K) = 
T>Gal(L x /K), where C x := = A*. 

By step 3, there exists a polynomial of degree n satisfying the relation (|A.1[) : 



A++) 


d 2 r P n - 


Y2 


(r + d r P n - 2nP n 


= o, 


A g A_ 


A+-) 


df.P n - 


h2 


(r - £) d r P n - 2nP n 


= o, 


Ae A_ 


A-+) 


d 2 r P n - 


h2 


{-r + *±!) d r P n + 2nP n 


= o, 


Ag A+ 


A-) 


d?P n - 


h2 


{-r - £) d r P n + 2nP n 


= o, 


Ae A. 



These polynomials exist for all A G A when their degrees are n G 2Z, while for 

n G 2Z + 1, they exist only for the cases A |_) and A ) with special conditions. 

Hence we have obtained the algebraic spectrum A = 2Z, where A++ = 4Z_ , A^ = 

2Z_, A_ + = 4Z+, A__ = 2Z. 

The possibilities for eigenfunctions, considering only A G 4Z, are given by 



A++) 


*n(r) 


= ^ +1 P 2n (r)e4, 


A G A_ 


A+-) 


*»(r) 


= f~ f2n(? , )e~, 


A G A_ 


A-+) 


*„(r) 


= r £+1 P 2 „(r)e^, 


A G A+ 


A-) 


*n(r) 


-r 2 


A G A. 



To obtain the point spectrum, we look for \& n satisfying the bound state conditions 

which is in only true for A G A |_. With the change r i— > y/^r, the point spectrum 

and ground state of the Schrodingcr equation with the 3D-harmonic oscillator po- 
tential are respectively Spec p (i/) = {E n : n G Z + }, for E n = 2nuj, where uj is the 
angular velocity, and 

The bound state wave functions are obtained as = p2 n '&o- Now, we can see 
that DGal(L /K) = B and £(H) = Vect(l). Since = P 2 n*o, for all A G A we 
have DG&\(L\/ K) = B and £(H - A) = Vect(l). In particular, BGa\(L x /K) = B 
and £(H - A) = Vect(l) for all A G Spec p (iJ), where A = %E. 

We remark that the Schrodinger equation with the 3D-harmonic oscillator po- 
tential, through the changes r i-> \ojr 2 and \t i— >■ y / r\E', fall in the class of Whittaker 
differential equations in which the parameters are given by 

(21 + 2,)uj + 2E 1 1 

K=- , U =-£+-. 

4a; ' P 2 4 



By Martinet-Ramis theorem, we have integrability when ±k ± /i is a half integer. 
These conditions coincides with our four sets A±±. 
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Coulomb potential: 

cf_ £(£+1) 
r + r 2 ' 4(1+1) 

The Schrodinger equation in this case can be written as 



dl^ = V „ ' h TTT 7T77 - £M 



V r 2 r 4(£+l) 2 
By the change r i— > 2 ^^" 1 ^ r we obtain the Schrodinger equation 

3 2 *=f^-^±lUl-AW A = l^£ 



and in order to apply Kovacic algorithm, we set 

r 2 r 

First we analyze the case for A = 1: we can see that this equation only could 
fall in case 2 or in case 4 of Kovacic's algorithm. We start discarding the case 
2 because by step 1 we should have conditions C2 and 003. So we should have 
E c = {2,4 + 41, -41} and = {1} and by step 2, we would have n g Z. Thus, 
n ^ Z + and D = 0, so the differential Galois group of this Schrodinger equation for 
A = 1 is SL(2,C). ' 

Now, we analyze the case for A ^ 1: we can see that this equation could fall 
in case 1, in case 2 or in case 4. We start discarding the case 2 because by step 1 
we should have conditions C2 and 003, so that we should have E c = {2, 4 + 4£, —4£} 
and Eoa = {0}. By step 2, we should have n = 2^GZ + ,so that D = {2£} and the 
rational function is 9 = but we discard this case because there could only 
exist one polynomial of degree 2£ for a fixed £, and for instance, there could only 
exist one eigenstate and one eigenfunction for the Schrodinger equation. 

Now, we only work with case 1, by step 1, conditions Ci and 003 are satisfied. 
Thus, 



1 



vT^A 

By step 2 we have the following possibilities for n e Z + and for A G A: 
A++) n = a+-a+ = -(£ + l)(l + 7 ^ % ), A = 1 - 

A+-) n = a+ -o^ = -^-+£, A = 1 - 



A-+) n = a~-a+ = (£ + !) bfb ~ 1 ), A = 1 



1+1 



2 



A__) n = a" -«o =^3.+^ A = l- |±i 
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We can see that A € A_ when A < 0, while A <G A+ when < A < 1. Furthermore: 



£<-l, Ae 



A+_) £>0, 



A € 



A_+) /eZ, Ae 



A__) £>0, 



A € 



A 

A_, 


1} ^ —n—2 
c — 2 




2 - 1 - 


A_, 


t — 2 


A+, 


< £< a=I 


A_, 


£<^ 


A+, 


£ > -1 


A_, 


£>^- 


A+, 


0<£<^- 



So the possible algebraic spectrum can be A = A ++ U A_| U A |_ U A , that is 



(5.2) A = <h - 



i+1 



neZ + U 1- 



l+l + n 

Now, for A £ A, the rational function uj is given by 

AeA++, R n 



£ + 1 
t-n 



n £ Z_i 



A++) 
A+-) 
A-+) 
A__) 



i+l+n ' r ' 



t(j+i) _ 2(^+1) _|_ ( e+i 



e+i _ I 

£ — n r 



Rn — 



l+l+n 



-n r 



x e A+ 

+ *±!, AeA_+, Rn 
AG A — , Rn 





7' 


l(t+l) 
r l 


2(1+1) 
r 


1(1+1) 

- 


2(£+l) 

7' 


_ 1(1+1) 


2(e+i) 

7' 



f+l + n J ' 



+ 



+ 



l+l + n I ' 



where R n is the coefficient of the differential equation = Rn^, which is inte- 
grablc for every n. 

By step 3, there exists a polynomial of degree n satisfying the relation (IA.1I) . 
A ++ ) ^»+ 2(^ + ^)^ + 2^(1 + ^)^ = 0, 



0. 



A + -) %Pn + 2 (ft - i) d r P n + « (l - P„ 



2(1+1) 



1+1 _ i 

i—n r 



8 P 



2(1+1) 



1 



i+1 



Pn 



0. 



These polynomials exist for every A £ A when neZ, but Pq = 1 is satisfied 

only for A G A |_. Hence we have confirmed that the algebraic spectrum A is given 

by equation (|5.2[) . 
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A. 



The possibilities for the eigenfunctions are given by 
-+) ¥»(r) = r e+1 P n (r)f n (r)e r , f n (r) = e***& , AG 



A+_) *„(r) = r"'P„(r)/„(r)e r , fn(r) = e~ 



X G 



A-+) * n (r) = r e+1 P n (r)f n (r)e- r , f n (r) = e?+^ , A G 



A__) *„(r) = r- l P n (r)f n (r)e- r , /»(r) = e 



S3' 



A G 



A 

i V — j 


(1 s -n-2 
t - 2 


A+, 


< i < 


A-, 


* — 2 


A+, 


< £ < ZL r r i 


A-, 


* ^- —n — 2 
fc — 2 


A+, 


f > -1 


A-, 


* — 2 


A+, 


< I < V" 



but ty n should satisfy the bound state conditions which is only true for A G A |_ n 

A + , so that we choose A |_ n A + = spec p (-ff), that is 



spec p (i?) = I 1 - 



1 



£ + n+ 1 



By the change r i— > 2 (i+i) r ' ^ ne P om t spectrum and ground state of the Schrodingcr 
equation with Coulomb potential are respectively 

Q 4 



Spec J H) = {E n :neZ+}, E n = - 



1 



1 



4 V^+i) 2 (^ + l + n) : 



and 



t+i 



e 2(*+i)' 



2(e+l) 

The eigenstates are given by ^/ ra = P n fn^0: where 

/„(r) = e 2(f + l + n )(f+l) . 

Now, we can see that DGal(L /A') = 1 and £ (H) = Vcct(l). Since \I>„ = P 2 «/n*o, 
for all A e A we have DGal(L A /A') = B and £(H - A) = Vect(l). In partic- 
ular, DGal(L A /AT) = B and £(H - A) = Vcct(l) for all E G Spec p (H), where 

4 

P C \ 

^ — 4(^+1)2 ^' 

We remark that, as in the three dimensional harmonic oscillator, the Schrodingcr 
equation with the Coulomb potential, through the change 



r i — y 



y/-4 {£+l) 2 E + c 4 



1 



falls in the class of Whittaker differential equation with the parameters given by 
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Applying Martinet- Ramis theorem, we can impose ±k±^ half integer, to coincides 
with our four sets A±±. 

Remark 5.5. By direct application of Kovacic's algorithm we have: 

• The Schrodingcr equation (|3.1j) with potential 

V(x) = ax 2 + -r 
x z 

is integrable for A £ A when 
-a = 0, b = n{p + l), ^eC, A = C, 
-o=l, 6 = 0, A G A = 2Z+ 1, 
-o=l, b = £(£+!), £ G Z*, A = 2Z + 1. 

• The only rational potentials (up to transformations) in which the elements 
of the algebraic spectrum are placed at the same distance, belongs to the 
family of potentials given by 

2 

V{x)= a k xk > °2^0. 

k— — 00 

In particular, the set A for the harmonic oscillator (a = 1, b = 0) and 3D 
harmonic oscillator (a = 1, b = £(£ + 1)) satisfies this. 

Proposition 5.6. Let C\ be the Schrddinger equation (|3.1[) with K = C(x) and 
Picard-Vessiot extension L\. If DGal(£o/-ftT) is finite primitive, thenDGal(L\/K) 
is not finite primitive for all A G A \ {0}. 

Proof. Pick V G C(x) such that £q falls in case 3 of Kovacic algorithm, 
then ouoo > 2. Assume t, s G C[x] such that V = |, then deg(t) > deg(s) + 2 
and V — A = a ~ . Now, for A ^ we have deg(.s — Xt) = dcg(t) and therefore 
o(V — A)oo = 0. So, for A 7^ 0, the equation C\ does not falls in case 3 of Kovacic 
algorithm and therefore DGal(LA/-?0 is not finite primitive. □ 

Corollary 5.7. Let C\ be the Schrddinger equation (13.11) with K = C(x) and 
Picard-Vessiot extension L\. //Card(A) > 1, then there is either zero or one value 
of A such that DGal(L^/AT) is a finite primitive group. 

Proof. Assume that Card(A) > 1. Thus, by proposition 15.61 the Schrodinger 
equation does not falls in case 3 of Kovacic's algorithm. □ 

We note that the study of the differential Galois groups of the Schrodingcr 
equation with the Coulomb potential had also been analyzed by Jean-Pierre Ramis 
( |51j ) using his summability theory in the eighties of the past century. 

6. Transcendental Potentials and the Algebrization Method 

In supersymmetric quantum mechanics, there exists potentials which are not 
rational functions and, for this reason, it is difficult to apply our Galoisian ap- 
proach from section [5] In this section, we address this problem and present some 
results concerning differential equations with non-rational coefficients. For these 
differential equations it is useful, whenever possible, to replace it by a new differen- 
tial equation over the Riemann sphere P 1 (that is, with rational coefficients). We 
will propose techniques to find such a change of variable and demonstrate, through 
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numerous examples, how our approach allows us to handle many examples from 
quantum mechanics. The equation over P 1 is called an algebraic form or algebriza- 
tion of the original equation. 

This algebraic form is already studied in the 19th century litterature (Liou- 
ville, Darboux), but the problem of obtaining an algebraic form (if it exists) of a 
given differential equation is in general not an easy task. Here we develop a new 
methodology using the concept of Hamiltonian change of variables. Such a change 
of variables allows us to compute the algebraic form of a large number of differen- 
tial equations of different types. In particular, for second order linear differential 
equations, we can apply Kovacic's algorithm over the algebraic form to solve the 
original equation. 

The following definition can be found in [9], 1691 170] ■ 

Definition 6.1 (Pullbacks of differential equations). Let £j £ ifi[<9 z ] and £2 £ 
be differential operators. The expression £2 <8> (d x + v) refers to the operator 
whose solutions are the solutions of £2 multiplied by the solution e~ J vdx of d x + v. 

• £2 is a proper pullback of £1 by means of / £ if 2 if the change of variable 
z = f{x) changes £1 into £2. 

• £2 is a pullback (also known as weak pullback) of £1 by means of / £ if 2 
if there exists v £ if 2 such that £2 <£> (d x + v) is a proper pullback of £1 
by means of /. 

In case of compact Riemann surfaces, the geometric mechanism behind the 
algebrization is a ramified covering of compact Riemann surfaces, see |44L 143] , 

6.1. Second Order Linear Differential Equations. Some results presented 
in this subsection also can be found in [4] §2 ]. 

Proposition 6.2 (Change of the independent variable, [4]). Let us consider the 
following eguation, with coefficients in C(z): 



and<C(z) >• L the corresponding Picard-Vessiot extension. Let (K,8) be a differen- 
tial field with C as field of constants. Let 9 £ K be a non-constant element. Then, 
by the change of variable z = 6(x), equation (|6.1| is transformed in 



Let Kq C K be the smallest differential field containing and C. Then equation 
(|6.2[) is a differential equation with coefficients in Kq. Let Kq Lq be the corre- 
sponding Picard-Vessiot extension. Assume that 

C{z) -> Kq, Z^9 

is an algebraic extension, then 

(DGal(Lo/A' ))° = (DGal(L/C(z)))°. 

Proposition 6.3. Assume that C x and C z are as in vrovosition 1 6. 2i Let tp be the 

transformation given by 



(6.1) 



Cz ■■= dly + a(z)d z y + b{z)y = 0, 




r — y o 0. 
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Then DG&l{L /K ) ^ DGal(L/K n L) C DGal(X/C(z)). Furthermore, if Kq n L 
zs algebraic overC(z), then (DGal{L / I<o))° ~ (DGal(L/C(ar)))°. 

Proof. By Proposition 16.21 the transformation (/? leads us to 

C(z) ~ p(C(«)) A' . 

We identify C(z) with <^(C(z)), so that we can view C(z) as a subfield of Kq and 
then by the standard Galois theory diagram (see [331 171] ). 

Lo . 

DGal(L /X ) 



DGal(L/(Ln*"o)) 




Lf]K 



C(z) 

we have 

DGal(L /K Q ) ~ DGal(L/i^ n L) C DGal(i/C(z)). 
If -FiTo H L is algebraic over C(z), then 

(DGal(Lo/#o))° - (DGal(L/C(z)))°. 

□ 

For the remainder of this section we write z = z(x) instead of 6. 

Remark 6.4 (Hard Algcbrization) . The proper pullback from equation (|6.2p to 
equation (|6.1|) is an algcbrization process. Therefore, we can try to algebrize any 
second order linear differential equations with non-rational coefficients (proper pull- 
back) if we can put it in the form of equation (|6.2[) . To do this, which will be called 
hard algebrization, we use the following steps. 

Step 1: Find a term of the form (d x z) 2 in the coefficient of y to obtain a 

candidate for d x z and then z. 
Step 2: Divide the coefficient of y by (d x z) 2 to obtain a candidate b(z) and 

check whether b G C(z). 
Step 3: Add (d x z)/d x z to the coefficient of d x y and divide the result by 

d x z to obtain a candidate a(z) and check whether a £ C(z). 



Let us illustrate this method of hard algebrization on the following example. 



Example. In }58[ p. 256], Singer presents the second order linear differential 
equation 

d l r ~ n 1 , iA r - ( lnx + !) 2r = °" 
x(lna; + l) 
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To algebrize this differential equation we choose (d x z) 2 = (In 2;+ l) 2 , so that d x z 
In x + 1 and for instance 



z = J (lna; + l)dx = xhix, b(z) = — 1. 
Now we find a(z) in the expression 
a(z)(\nx + 1) - 



x(hix + l) x(lnx + l)' 

obtaining a(z) = 0. So that the new differential equation is given by d 2 y — y = 0, 
in which y(z(x)) = r(x) and one basis of solutions of this differential equation is 
given by (e z ,e~ z ). Thus, the respective basis of solutions of the first differential 
equation is given by (e xlnx , e _x ). 

In general, this method is strongly heuristic because the quest of z — z(x) 
in b(z)(d x z) can be purely a lottery, or simply there may not exist z such that 
a{z),b{z) € C(z). For example, the equations presented by Singer in [581 p. 257, 
261, 270] and given by 

„ 9 =p2xln 2 x =F 2x\nx — 1 „ — 2xln 2 x — 3xlna; — x =F 1 

d 2 x r + ± ^_ d x r + = 0, 

3 2 H ^r— r = 0, (x 2 In 2 x)<9 2 r+(a; In 2 x— 3x\nx)d x r+3r = 0, 

4x In x 4x In a; 

cannot be algebrized systematically with this method, although they correspond to 

pullbacks (not proper pullbacks) of differential equations with constant coefficients. 

In [14] , Bronstein and Fredet developed and implemented an algorithm to solve 
differential equations over C(x, e$ f) without algebrizing them, see also |26j . As an 
application of proposition 16.21 we have the following resultQ- 

Proposition 6.5 (Linear differential equation over C(x,ef f ), H). Let f e C(x) 
be a rational function. Then, the differential equation 



r = 0, 



(6.3) d 2 x r -(f+f*d- fef 'a (e' >) W + (/ (e' ')) * b (V ') 

is algebrizable by the change z = e$ ? and its algebraic form is given by 
d 2 y + a(z)d z y + b(z)y = 0, r(x) = y(z(x)). 
Proof. Assume that r(x) = y(z(x)), and z = z(x) = e^ dx . We can see that 

— *-*(S)-^(*-'-(¥))- 

Replacing in d 2 y + a(z)d z y + b(z)y = we obtain equation (|6 .3[) . □ 



This result is given in [4] §2], but we include here the proof for completeness. 
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Example. The differential equation 

figr - (x + ~ - 2xe x2 j d x r + Xx 2 e x \ = 0, 

2 

is algebrizable by the change z = e~ and its algebraic form is given by 

d 2 z y + 2zd z y + Xy = 0. 

Remark 6.6. According to proposition ^. 5[ we have the following cases. 

(1) / = n^f^, for a rational function h, n € Z+, we have the trivial case, both 
equations are over the Riemann sphere and they have the same differential 
field, so that does not need to be algebrized. 

(2) / = for a rational function h, n € Z + , (|6.3|) is defined over an 
algebraic extension of C(x) and so that this equation is not necessarily 
over the Riemann sphere. 

(3) / 7^ f° r an y rational function h, q e Q, (|6.3p is defined over a 
transcendental extension of C(x) and so that this equation is not over the 
Riemann sphere. 

To algcbrizc second order linear differential equations may be easier when the 
term in d x r is absent and the change of variable is Hamiltonian. 

Definition 6.7 (Hamiltonian change of variable, [4]). A change of variable z = z(x) 
is called Hamiltonian if [z(x), d x z(x)) is a solution curve of the autonomous classical 
one degree of freedom Hamiltonian system 

/ d x z = d w H . w 1 

1 d x w = -d z H Wlth H = H ^ W) = T + V{zl 

for some V € K. 

Remark 6.8. Assume that we algebrize equation (|6.2[) through a Hamiltonian 
change of variables z = z(x), i.e. V G C(z). Then A' = C(z, d x z, . . .), but we have 
the algebraic relation 

{d x zf = 2h-2V(z), h = H(z, d x z) e C, 

so that Kq = C(z,d x z) is an algebraic extension of C(z). By proposition 16. 2\ the 
identity connected component of the differential Galois group is preserved. On 
the other hand, we can identify a Hamiltonian change of variable z = z(x) when 
there exists a £ K such that (d x z) 2 = ot(z). Thus, we introduce the Hamilton- 
ian algebrization, which corresponds to the algebrization process done through a 
Hamiltonian change of variable. 

The following result, which also can be found in [4J §2], is an example of 
Hamiltonian algebrization and correspond to the case of reduced second order linear 
differential equations. 

Proposition 6.9 (Hamiltonian Algebrization, [4j). The differential equation 

d 2 r = q(x)r 

is algebrizable through a Hamiltonian change of variable z = z{x) if and only if 
there exist /, a such that 

— , - e C(z), where f(z(x)) = q(x), a(z) = 2(H - Viz)) = (d x z) 2 . 
a a 



GALOISIAN APPROACH TO INTEGR ABILITY OF SCHRODINGER EQUATION 35 



Furthermore, the algebraic form of the equation d 2 r = q{x)r is 

(6.4) d 2 y+^-—d z y~-y = 0, r(x) = y(z(x)). 

2 a a 

Remark 6.10 (Using the Algebrization Method). The goal is to algebrize the 
differential equation d 2 r = q{x)r. We propose the following steps. 

Step 1: Look for a Hamiltonian change of variable z = z{x) and two func- 
tions / and a such that q(x) = f(z(x)) and (d x z(x)) 2 = a(z(x)). 
Step 2: Verify whether or not f(z)/a(z) £ C(z) and d z a{z)/ a{z) £ C(z) to 

see if the equation d 2 r = q(x)r is algebrizable. 
Step 3: If the equation d 2 r = q(x)r is algebrizable, its algebrization is 

dzV + l — d z y - -y = 0, y(z(x)) =r(x). 
2 a a 

When we have algebrized the differential equation d 2 r = q(x)r, we study its inte- 
grability, its eigenring and its differential Galois group. 

Examples. Consider the following examples. 

• Given the differential equation d 2 r = /(tanx)r with / £ C(tanx), we can 
choose z = z(x) = tana; to obtain a(z) = (1 + z 2 ) 2 , so that z = z(x) 
is a Hamiltonian change of variable. We can see that — £ C(z) and 
the algebraic form of the differential equation d 2 r = /(tanx)r with this 
Hamiltonian change of variable is 

2z f(z) 

i + 2^ (l + z i y 

• Given the differential equation 



o2 VI + x 2 + x 2 

o x r = r, 

x 1 + x 2 



we can choose z = z(x) — yl + x 2 to obtain 

f(z) = - 2 , a{z) = 

so that z = z(x) is a, Hamiltonian change of variable. We can see that 
^ £ C(z) and the algebraic form for this case is 

2 



^ z -\- z 1 

^ + —(—2 — T\ d *y 2 — r~ V = °> y(Vl + a?)=r(x). 

z{z z — 1) z- 2 — 1 

We remark that in general the method of Hamiltonian algebrization is not an 
algorithm, because it may not be obvious to obtain a suitable Hamiltonian H satis- 
fying definition 16. 71 We present now a particular case of Hamiltonian algebrization 
considered as an algor ithnfl 

Proposition 6.11 (Hamiltonian Algebrization Algorithm, [4]). Consider q(x) = 
g(zi, ■ ■ ■ ,z n ), where Zi = e XiX , Xi £ C*. The equation d 2 r = q(x)r is algebrizable 
if and only if. 

^eQ*, l<i<n, l<j<n, g£C(z). 



•^Proposition [6TTT1 is a slight improvement of a similar result given in [4] §2], Furthermore, 
we include the proof here for completeness. 
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Furthermore, A; = c,A, where A £ C* and Cj € Q* and /or the Hamiltonian change 
of variable 

z = e~ , where a = — , pi,q, e Z*, gcd (p*, <&) = 1 and o = lcm(gi, ■ ■ • , g„), 
t/ie algebrization of the differential equation d^r = q(x)r is 

o z y+-o z y-q y = 0, rm = — , = r(x). 



Proof. Assuming Aj/Aj = Cy € Q* we can see that there exists A G C* and 
Cj S Q* such that Ai = Ac;, so that 



Now, setting z = z{x) = ei x we can see that 

f(z) = g(z m \...,z m -), mi - 



2 ^2 



A 2 z 



As q\qi, we have mj € Z, so that 



, - e C(z) 

a a 



and the algebraic form is given by 

dlv+-d zV -q 2 ^ ^ >-y = 0, y(z(x)) = r(x). 

Z A Z Z Z 



□ 



Remark 6.12. Propositions 16.91 and 16.111 allow the algebrization of a large number 
of second order differential equations, see for example |49j . In particular, under 
the assumptions of proposition 16. Ill we can algcbrize automatically differential 
equations with trigonometrical or hyperbolic coefficients. 

Examples. Consider the following examples. 
• Given the differential equation 

Q , e** + 3e-§* -2ei x 1 2 x 5 x , , 3 

o x r = — r r, Ai = -, A 2 = --, A 3 = -, A 4 = 1, A 5 = --, 

e -\- e 2 ^ o4 z 

we see that Xi/Xj € Q, A = 1, q = lcm(l, 2, 3, 4) = 12 and the Hamiltonian 

change of variable for this case is z = z(x) = e~^ x . We can see that 



a ( z ) = 771 z ' = 1 2 , _is » ,-eC(z) 

144 z L/S + z 18 a a 

and the algebraic form is given by 

d 2 z y+ -d z y- 144 Z + ^ 2z y = 0, y(e*f a ) = r(x). 
• Given the differential equation 

figr = (e 2V5x + e"^ - e 3a >, Ai = 2^2, A 2 = -V2, A 3 = 3, 

we see that A1/A2 G Q, but A1/A3 ^ Q, so that this differential equation 
cannot be algebrized. 
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We remark that it is also possible to use the algebrization method to transform 
differential equations when either the starting equation has rational coefficients or 
the transformed equation does not have rational coefficients. 

Examples. As an illustration, we treat the following examples. 

• Consider the following differential equation: 

2 x A + 3x 2 — 5 

cLr = - r = 0. 

x 2 + 1 

We can choose z = z(x) = x 2 so that a — 4z and the new differential 
equation is 

d 2 y + —d z y -— ^ -r— V = 0, y(x 2 ) = r(x). 

zy 2z y 4z(z + l) y ' yy ' w 

• Consider the Mathieu differential equation d 2 r = (a + bcos(x))r. We can 
choose z(x) = ln(cos(x)) so that a = e~ 2z — 1 and the new differential 
equation is 

I ae 2z _|_ 5g3z 

d *y - 1 — 2i d *y — i — = °> y(Mc°s(x))) = r ( x )- 

1 — e^ z 1 — e 

Recently, the Hamiltonian algebrization (propositions 16.91 and I6.11[) has been 
applied in [H [2l |4] to obtain non-intcgrability results. 

6.2. The Operator d z and the Hamiltonian Algebrization. The gener- 
alization of proposition 16 . 21 to higher order linear differential equations is difficult. 
But, it is possible to obtain generalizations of proposition 16.91 by means of Hamil- 
tonian change of variable. We recall that z = z(x) is a Hamiltonian change of 
variable if there exists a such that (d x z) 2 = a(z). More specifically, if z = z(x) is 
a Hamiltonian change of variable, we can write d x z = y/a, which leads us to the 
following notation: d z = yfad z . 

We can see that d z is a derivation because satisfy d z (f + g) = d z f + d z g and 
the Leibnitz rules 

S.(f. g ) = 8.f. g + f-8,g, 8.U) = 'fa . 

\9J 9 

We can notice that the chain rule is given by d z (fog) = d g f ogd z (g) ^ d g f ogd z (g). 
The iteration of d z is given by 

fi£ = l, d z = V^d z , = V^dJ?- 1 = V^d z (. . . (V^d z )) . 

S v ' 

n times */ad z 

We say that a Hamiltonian change of variable is rational when the potential V G 
C(x) and for instance a € C(x). For the remainder of this paper, we let d z = \fad z 
where z = z(x) is a Hamiltonian change of variable and d x z = \fa. In particular, 
d z = d z = d x if and only if •Jot = I, i.e., z = x. 

Theorem 6.13. Consider the systems of linear differential equations [A] and [A] 
given respectively by 

8 X Y = -AY, d z Y = —AY, A=[a ij ], A = [ay], Y = [yn], Y=[y a ], 
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where aij e K = C(z(x),d x (z(x))), ay <G C(z) C A = C(z,^/a), 1 < i < n, 
1 < j ' < n j a ij( x ) = ®ij( z { x )) an d Un(x) = yn(z(x)). Suppose that L and L are 
the Picard-Vessiot extensions of [A] and [A] respectively. If the transformation tp 
is given by 

yn(x) h> ya(z(x)) ' 
d x i-» d z 

then the following statements hold. 

. K~K, (K,d x ) ~ (K,d z ). 

• DGal(L/A) ~ DGal(L/A) C DGal(Z/C(z)). 

• (DGal(L/A))° ~ (DGal(L/C(z))°. 

• £{A)~£{A). 

Proof. We proceed as in the proof of proposition ^. 31 As z = z(x) is a rational 
Hamiltonian change of variable, the transformation if leads us to 

C(z) ~ ^(C(z)) A, A ~ A, C(z) A, (A, ft.) ~ (A, §,) 

that is, we identify C(z) with ip(C(z)), and so that we can view C(z) as a subfield 
of A" and then, by the standard Galois theory diagram (see [331 171| ). 

L . 



DGal(L/_R") 



DGal(L/if) 




C(z) 

so that we have DGal(L/A) ~ DGal(i/A) C DGal(L/C(z)), (DGal(i/A))° ~ 
(DGal(L/C(z)))°, and £{A) ~ 8(A). □ 

We remark that the transformation <p, given in theorem l6.131 is virtually strong 
isogaloisian when y/a ^ C(z) and for y/a £ C(z), if is strong isogaloisian. Further- 
more, by cyclic vector method (see [65]), we can write the systems [A] and [A] in 
terms of the differential equations C and C. Thus, C is the proper pullback of C 
and £(£) ~ £ (£). 

Example. Consider the system 

r a _ 2\/2 
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Through the Hamiltonian change of variable z = e x , and with yfa = z, it is 
transformed into the system 



[A] := 



'^(i-z 2 y 



f) Z 2^2 z z 2 -l 

One solution of the system [^4] is given by 

/V2 

1 / ~ 

z 2 + 1 

and the corresponding solution for the system [A] is 

1 / 2 



2; 

—z 



e 



Remark 6.14. The algebrization given in proposition 16.91 is an example of how 
the introduction of the new derivative d z simplifies the proofs and computations. 
Such proposition is naturally extended to d x y + ad x y + by = 0, using ip to obtain 
d 2 y + ad z y + by = 0, which is equivalent to 

(6.5) ad 2 y + + V^a) d z y + by = 0, 



2 

where y(x) = y(z(x)), a{z{x)) = a(x) and b(z(x)) = b(x). 

In general, for y(x) = y(z(x)), the equation F(d x y, . . . , y, x) = with co- 
efficients given by a,i k (x) is transformed into the equation F(d"y, . . . , y, z) = 
with coefficients given by ai k {z), where a,i k {x) = ai k (z(x)). In particular, for 
y/a, a,i k G C(z), the equation F(d z y 1 . . . , y, z) = is the Hamiltonian algebriza- 
tion of F(d x y, . . . , y, x) =0. Now, if each derivation d x has even order, then a and 
a,i k can be rational functions when we algebrize the equation F{d x y, . . . , y, x) = 0, 
where ai k £ C(z(x), d x z{x)). For example, that happens for linear differential 
equations given by 

d 2 x n y + a n -x{x)d 2 x n ~ 2 y + ... + a 2 {x)d A x y + ai(x)d*y + a (x)y = 0. 

Finally, the algebrization algorithm given in proposition 16 . 1 ll can be naturally ex- 
tended to a differential equation of the form 

which, through the change of variable z = e^ x , is transformed into 

F{d:y 1 d n z - 1 y,...,d z y,y,z) = Q- 

In particular, we may consider the algebrization of Riccati equations, higher order 
linear differential equations and systems. 
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Examples. The following corresponds to some examples of algebrizations for dif- 
ferential equations given in 1581 p. 258, 266]. 

(1) The equation C := d 2 y + (~2e x — l)d x y + e 2x y = with the Hamiltonian 
change of variable z = e x , \fa = z,a = —2z — 1 and b = z 2 is transformed 
into the equation C := d 2 y—2d z y+y = which can be easily solved. Bases 
of solutions for C and C are respectively given by (e z , ze z ) and (e 6 * , e I e e "). 
Furthermore A = C(e x ), K = C(z), L and L are the Picard-Vcssiot 
extensions of C and C respectively. Thus, DGal(L/A") = DGal(L/A'). 

(2) The differential equation 

n2 -24e a; -25 n 206* 

£: =^ + ^T^^ + 4^T5 y = 

with the Hamiltonian change of variable z = e x , \fa = z, 

_ -24z-25 - 20z 
and o = 



4z + 5 4z + 5 

is transformed into the equation 

? ^ -20(z + 1) ^ 20 ^ 
£: =^ + x(4z + 5) ^ + ^Tsf = °' 

which can be solved with Kovacic algorithm. A basis of solutions for C 
is (z + l,z 5 }, so that a basis for C is (e x + l,e 5x ). Furthermore K = 
C(e x ), K = C(z), L and L are the Picard-Vessiot extensions of C and C 
respectively. Thus, DGal(L/A) = DGal(L/A) = e. 

Remark 6.15 (Algcbrization of the Riccati equation). The Riccati equation 

(6.6) d x v = a(x) + b(x)v + c(x)v 2 . 

Through the Hamiltonian change of variable z = z(x), it is transformed into the 
Riccati equation 

(6.7) d z v = -^(a(z) + b(z)d + 'B(z)d 2 ), 



where v(x) = v(z(x)), a(z(x)) = a(x), b(z(x)) = b(x), c(z(x)) = c(x) and \J a(z(x)) = 
d x z(x). Furthermore, if y/a, a, b, c G C(x), equation (|6.7[) is an algcbrization of 
equation (|6.6|) . 

Example. Consider the Riccati differential equation 

C := d x v = ( tanh x : — J v + (3 tanh x — 3 tanh 3 x) v 2 , 

\ tanhxy v ' 

which through the Hamiltonian change of variable z = tanhx, for instance ^fa = 
1 — z 2 , is transformed into the Riccati differential equation 

C := d z v = — v + 3zv 2 . 

z 

One solution for the equation C is 

; r-, where c is a constant, 

2(3,2 — c) 
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so that the corresponding solution for equation C is 

1 



tanh ir(3 tanh x — c) 

The following result is the algcbrizcd version of the relationship between the 
eigenring s of systems and operators. 

Proposition 6.16. Consider the differential fields K, K and consider the systems 
[A] and [A] given by 

d x X = -AX, d z X = -AX, d z = s/ad z , A = [ay], A = [%], a tj G K, a tJ G K, 

where z = z(x), X(x) = ~K.(z(x)), aij(z(x)) = a i j{x), then their eigenrings satisfy 
£{A) ~ £{A). In particular, if we consider the linear differential equations 

n— 1 n—1 

C:=d^y + J2 a hd^y = and £:=^|^a^!/ = 0, 

fc=0 fc=0 

where z = z(x), y(x) = z{{x)), a.k(z(x)) = ak(x), G i^, G If, i/ien £ (£) ~ 
£(£), where C := £y = and £ := £j/ = 0. Furthermore, assuming 

fan •■• Pi»\ / 1 ... 

La- 

\Pnl ■■■ Pnn) \—OQ ~a\ ... —On-lj 

then 

£(£) = j&i*^ -1 ■ d x P = PA- AP, Pik G if j , 

zf and onh/ if 

£(£) = jx>lft^ _1 : d z P = PA-AP, p lk G £ j . 

Proof. By theorem iH we have K ~ £(A) ~ £(A) and £(£) ~ £(£). 
Using the derivation d z and by induction on lemma [6~51 we complete the proof. 

□ 

Examples. We consider two different examples to illustrate the previous proposi- 
tion. 

• Consider the differential equation C\ := d x y — (1 + cos a; — cos 2 x)y = 0. 
By means of the Hamiltonian change of variable z = z(x) = cos a;, with 
\Ja = —yl — z 2 , L\ is transformed into the differential equation 

Z „ _ 1 + 2- 2 2 

1 — 2 d *y — i — 2 

1 — 2 Z 1 — Z z 



£i ■= 9 z y- ; - 2 ^y- - ' o~ ^ = °- 

>^ — — — 



Now, computing the eigenring of £i we have £(£\) = Vect(l), therefore 
the eigenring of £i is given £(£i) = Vect(l). 
• Now we consider the differential equation £ 2 := d x y = (e 2x + |) y. By 
means of the Hamiltonian change of variable z = e x , with *Ja = x, £2 is 
transformed into the differential equation 

£ 2 :=d 2 z y+- z d z y-(l + ^)y = 0. 
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Now, computing the eigenring of £2 we have 
£(&) = Vect (l,-2 ( *=±) d z - *^*) = Vect ( 1, -2 / 



therefore the eigenring of £2 is given by 

£(£2) = Vect (l, -2 (^—^] d x - 



The same result is obtained via matrix formalism, where 

with P G M<rf(2, C(e x )) and P G Mat(2,C(z)). 

6.3. Applications in Supersymmetric Quantum Mechanics. In this sub- 
section we apply the derivation d z to the Schrodingcr equation = X$!, where 
H = — d% + V(x), V G K. Assume that z = z(x) is a rational Hamiltonian change of 
variable for = A if, then K = C(z(x), 9 x z(o;)). Thus, the algebrized Schrddinger 
equation is written as 

(6.8) #§ = A§, P = -§* + V^», d 2 =ac\ 2 + i<9 z a<9 2 , A' = C(z,Va). 
The reduced algebrized Schrddinger equation is given by 

V(*)=V + XM 

(6.9) H$ = A$, H = a(z) (-<9 2 + V(z)) , V = c4W + W 2 , 

^ ~ 4 a( Z ) ' 

The eigenfunctions ^, $ and $ corresponding to the operators H, H and H are 
related respectively by 

$(z(x)) = ^fa^{z{x)) = ^a^{x). 

In order to apply the Kovacic's algorithm, we only consider the algebrized 
operator H, whilst the eigenrings will be computed on H. Also it is possible to 
apply the version of Kovacic's algorithm given in reference |62j to the algebraized 
operator H. The following results are obtained by applying Kovacic's algorithm to 
the reduced algebrized Schrodingcr equation (equation (|6.9[1 ) H$ = A3>. 

Proposition 6.17. Let L.\ be the Picard-Vessiot extension of the reduced algebrized 
Schrddinger equation H$ = A$ with a, V G C[z]. If dega < 2 + degV then, for 
every A 6 A, DGal(L/A") is not finite primitive. 

Proof. Suppose that dega = n and degV = m. The reduced algebrized 
Schrodingcr equation H<I > = A$ can be written in the form 

.^o ^ Aad'ia — 3(d z a) 2 + l6a(V - A) 

locr 

As to > n — 2 we have 0(7-00) = n — to < 2, which does not satisfy the condition 
(00) of case 3 of Kovacic's algorithm, therefore DGal(L/A') is not finite primitive 
for any A G A. 
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□ 

Proposition 6.18. Let Li\ be the Picard-Vessiot extension of the reduced algebrized 
Schrodinger equation H<f> = A$ with a € C[z], V € C(z). If o(V)oo < 2 — dega 
then, for every A G A, DGal(L>,/ K) is not finite primitive. 

Proof. Suppose that V = s/t, where s and t are coprime polynomials in 
C(z). Assume that dega = n, degs = m and degi = p. The reduced algebrized 
Schrodinger equation H$ = A$ can be written in the form 

„o ^ Atad 2 z a - 3t(d z a) 2 + lQa(s - Xt) 
a $ = r$ r = - 

z ' 16ta 2 

As m > n + p — 2, we have o(r oc ) = p + n — m < 2, which does not satisfy condition 
(oo) of case 3 of Kovacic's algorithm. Therefore DGal(L/A^) is not finite primitive 
for any A € A. 

□ 



Remark 6.19. In a natural way, we obtain the algebrized versions of Darboux 
transformation, i.e., the algebrized Darboux transformation, denoted by DT. The 
notation DT„ stands for the n-th iteration of DT; the notation CI n denotes the 
n-th iteration of the algebrized Crum iteration, where the algebrized wronskian is 
given by 

y i • • • Vn 



W(yi,...,y n ) 



yi 



d 1 : 



y, 



In the same way, we define algebrized shape invariant potentials, algebrized superpo- 
tential W , algebrized super symmetric Hamiltonians H±, algebrized supersymmetric 



partner potentials V±, algebrized ground state ^ 
>(-) 



(-) 



, algebrized wave 



functions ^ , algebrized raising and lowering operators A and A 1 . Thus, we can 
rewrite entirely section [3] using the derivation d z . 



The following theorem shows us the relationship between the algcbrization and 
the Darboux transformation. 



Theorem 6.20. Given the Schrodinger equation C\ := if_\&( _ J = X$f(', the 
relationship between the algebrization tp and Darboux transformations DT, DT with 
respect to C\ is given by DTyj = i^DT, that is T)Tip(C) = ipDT(C). In other words, 
the Darboux transformations DT and DT are intertwined by the algebrization tp. 

PROOF. Consider the equations C\ := = A* (_) , C\ := = 

A$(~), C x := = A*(+) and Z := H + ¥+) = A§W, where the Darboux 

transformations DT and DT are given by DT(£) = C, DT(£) = £, 

V- ^ v+ y_ ^ v+ 

DT : , DT : 
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and <p(£x) = £x, where the algebrization tp is given as in theorem 16. 131 Then the 
following diagram commutes 



£x *■ £ 



£ 



5t Cx 



so DT<^(£) = y>DT(£) i.e. £ = £. 



□ 



To illustrate DT we present the following examples. 



Examples. Consider the algebrized Schrodinger equation > = A\&' ' with: 

jij. Taking Ai = 1 and = Jf±f, 



\/a(z) = \J z 1 — 1 and 
we have DT(F_) = = ^ and 



DT^" 5 ) 



(+) 



= \/z 2 - 



(-) 

A 



1 



:VL7 



(") 



where ^\ is the general solution of ' = A^ ' for A 7^ 1. 

The original potential corresponding to this example is given by V_ (#) 
c ° sh x 1 and for Ai = 1 the particular solution is given by sl " h x , . 

cosna: — I x * J- J cosna; — 1 

™ shT _ and DT(*^ - 



Applying DT we have DT(VL) = V+(x) 



(+) 



9x* 



(-) 

A 



A 

• -v/ci = —z, V- (z) = z 
DT(t/L) = v: - - 2 



smn a: A 



cosh a:+l 



z. Taking Ai = and * 



+- z and DT(§^ ] ) 



we have 



(-) 

A ' 



where ^ x ' is the general solution of ' = A$' ' for A 7^ 0. This 

example corresponds to the Morse potential V-(x) = e~ 2x — e~ x . 

To illustrate the algebrized Crum iteration CI„, we present the following ex- 
ample, which is related with the Chebyshev polynomials. 



Example. Consider ^Ja. = — \/l — z 2 , V = with eigenvalues and eigenfunctions 
Ai = 1, A2 = 4, $1 = z, $4 = 2z 2 - 1, $„2 = T„(z), where T n (z) is the Chebyshev 
polynomial of first kind of degree n. The algebrized Wronskian for n = 2 is 



W(z, 2z 2 - 1) = - v/l - z 2 (2z 2 + 1), 
and by algebrized Crum iteration we obtain the potential 

CI 2 (tO = % = ((2z 2 - 1)<9 2 + zd z )lnl?(z,2z 2 - 1) 
and the algebrized wave functions 

W{z,2z 2 ~\,T n ) 



CI 2 ($a) = $f 



J¥(z,2z 2 - 1) 
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In a natural way, we introduce the notion of algebrized shape invariant po- 
tentials V n +i(z, a n ) = V n {z,a n +i) + R(a n ), where the energy levels for n > are 
given by E n = R{a\) + • • • R(a n ) and the algebrized eigenfunctions are given by 
^n(0'i) = A* (z, cti) • • ■ A^(z, a n )'^o(z, a n ). To illustrate the algebrized shape invari- 
ant potentials and the operators A and A\ we present the following example. 

Example. Let y/a = 1 — z 2 and consider the algebrized super potential W{z) = z. 
Following the method proposed in remark T4.101 step 1, we introduce fi € C to 
obtain W(z;n) = fiz, and 

V± {z; n) = W 2 {z- n) ± 8 z W{z; fx) = T l)z 2 ± fi. 

Thus V+(z; ao) = ao(ao — l)z 2 + ao and V-(z; ai) = a\{a\ + l)z 2 — oi. By step 2, 

d z {V+ (z;ao) - VL(«;oi)) = 2z(l - z 2 )(a (a - 1) - 01(01 + 1)). 

By step 3, we obtain 

ai(ai + 1) = a (a - 1), a 2 - a 2 Q = -(ai + a ) 

and, assuming a\ 7^ ±ao, we have a\ = /(ao) = ao — 1 and R{a\) = 2ao + 1 = 
(ao + l) 2 — ag = a 2 — «q. This means that the potentials V± are algebrized shape 
invariant potentials where the n-th energy level E = E n is easily obtained: 

n n 

E n = ^2 R ( a k) = X] ( a k ~ a k-l) = a n ~ a 0: a « = /™( a o) = a + 71. 
fc=l fc=l 

Now, the algebrized ground state wave function of V—(z,ao) is 

$0 = e J 



Finally, we can obtain the rest of eigenfunctions using the algebrized raising oper- 
ator: 

*„(z,a ) = A\z,a )A\z,a 1 ) ■ ■ ■ A f (z, a n -i)*o(z, a n ). 
This example corresponds to the Poschl- Teller potential. 

Now to illustrate the power of Kovacic's algorithm with the derivation d z , we 
study some Schrodinger equations for non-rational shape invariant potentials. 

Morse potential: V(x) — er 2x — e~ x . 

The Schrodinger equation H^f = A\I> is 

a 2 * = (e- 2:c - e~ x - A) 
By the Hamiltonian change of variable z = z(x) = e~ x , we obtain 

z 2 - z - i 

a{z) = z 2 , V(z) = z 2 -z, V(z) = 



z" 

Thus, K = C(z) and K = C(e x ). The algebrized Schrodinger equation = X^> 
is 

z 2 <9 2 § + zd^ -(z 2 -z- A)$ = 
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and the reduced algebrized Schrodinger equation H<J> = A$ is 

z 1 - z - \ - A 



cr$ = r$, r 

This equation could only fall under case 1, in case 2 or in case 4 (of Kovacic's 
algorithm). We start by analyzing the case 1: by conditions C2 and 003 we have 
that 

[Vr] = 0, c*o = 1±2 2 v/ ~^ , [y/rloo = 1 > and a ^ == F^- 

By step 2 we have the following possibilities for n 6 Z + and for A 6 A: 

,2 



,2 



A+ 


+) 




- a o 




A = 


-(n 


A+ 


-) 


« = a i 


- a 


= -1 + x/=A, 


A = 


-(n 


A_ 


+) 








A = 


-n 2 , 


A_ 


-) 




-«o 




A = 


-n 2 . 



We can see that A € A_ = {— rt 2 : n € Now, for A 6 A, the rational function 

lu is given by: 

A ++ ) w = i+3±fn A e A ++ , r„ = 4 " 2 ^"+ 3 + ^ + 1, 

A+ _) w = i_l±|n ; A e A + _, r n = 4 " 2 + 8 2 "+ 3 - + 1, 

A_+) w = -l + ^2, A e A_ + , r n = ^i-^i + l, 

A__) w = -l + i=fa, A e A — , r „ = + + L 
where r n is the coefficient of the differential equation 9 2< & = r„$. 

By step 3, there exists a polynomial of degree n satisfying the relation (|A.1[) . 
A ++ ) ff 2 z P n + 2 (1 + 3±3») a 2 p„ + ^±21 p n = 0, 

A+_) 9 z 2 P„ + 2(l-i±^)S,P n + ^Pn = 0, 
A_+) 5 2 P„ + 2(-l + ii|2)a 2 P„+ 2 i^P n = 0, 

A—) d 2 z P n + 2(-l + ^)d z P n + ^P n = 0. 

These polynomials only exist for n = A = 0, with A € A |_ U A So the 

solutions of = 0, = and H$ = are given by 

The wave function "tp satisfy the bound state conditions, which means that is 
ground state (see |23j ) and € Spec p (P). Furthermore, we have 

DGal(L /K) = DGal(L /£) = DGal(L /C(z)) = B, 
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£(H) = 8(H) = £(H) = Vect(l). 
We proceed with case 2. The conditions c 2 and 003 are satisfied. So we have 



(2, 2 + 4a/ = A, 2 - 4-\/ = A| and = {0}, 



and by step two. we have that 2 ± V — A = m G Z+, so that A = — f 12 ^ 1 ) and the 
rational function 8 has the following possibilities 

8 + = 2 -±^, 6.=-™. 

z z 

By step three, there exist a monic polynomial of degree m satisfying the recurrence 
relation (|A.2|) : 

fl ~1 £)3 p 1 3m+6 a2 p 4z 2 — 4z — 2m 2 — 7m — 6 a p 4mz+82-4m-6 p n 

C_|_J O z r m -\ j — O z r m O z f m "m — U, 

33p 3m o2 p 4z 2 — 4z — 2m 2 — m a p i 4mz— 4m — 2 p n 

1 7 — ; <J z Jr m - u z r m z 2 u z r rn T z i r rn — {J - 

We can see that, for m = 1, the polynomial exists only in the case 0_, with 
Pi = z — 1/2. In general, these polynomials could exist only in the case 8- with 
m = 2n — 1, ?7 > 1, that is A G {— n 2 : n > 1}. 

By case one and case two, we obtain A = {— n 2 : n > 0} = Spec p (i?). Now, 
the rational function <f> and the quadratic expression for uj are 



2ri - 1 d z P 2n -i 2 ~M ± \/M 2 - 4N 

z Pan-l 2 

where the coefficients M and iV are given by 

71/ = 2n ~ 1 _ d ^-i N = n2 - n +I _ -2 ^ d 2p 2ni ^ ^ 

Z -P2n-1 z2 Z ^2n-l 

Now, A = M 2 - AN ^ 0, which means that H$ = -n 2 $ with n e Z+ has two 
solutions given by Kovacic's algorithm: 

yfzP n e- z yfzP n -\e z 

*l.n = , $2,n = ~ ■ 

The solutions of H*S? = —n 2 ^ are given by 

Wl.n — — , W 2 ,„ — , 



and therefore, the solutions of the Schrodingcr equation H^f = — n 2 "^ are 

* i,» = P n e~ e ~*e nx , * 2>n =Pn-ie e ~*e nx , P n =P n oz. 
The wave functions = ty n satisfy the conditions of bound state and, for 
77 = 0, this solution coincides with the ground state presented above. Therefore we 
have 

$n = <f>ofnPn, $n = ^ofnPn, fn(z) = — . 

Thus, the bound states wave functions are obtained as 

*„ = *o/nP„ fn(x) = f„(e~ x ) = e nx . 
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The eigenring s and differential Galois groups for n > satisfies 

DGal(L„/A) = DGal(L„/A) = DGal(L„/C(z)) = G TO , 

dime £ (H + n 2 ) = dime £(H + n 2 ) = dim c £(H + n 2 ) = 2. 
We remark that the Schrodinger equation with Morse potential, under suitable 
changes of variables |38j , may be transformed to a Bessel differential equation. 



It is known that equations deriving from Eckart, Rosen-Morse, Scarf and Poschl- 
Tcllcr potentials may be mapped, under suitable transformations, to Hypcrgcomet- 
ric equations. These potentials are inter-related by point canonical coordinate 
transformations (see [161 p. 314] ), so that A = C because Poschl- Teller potential 
is obtained by means of Darboux transformations of V — I |42l 152] ). We consider 
some particular cases of Eckart, Scarf and Poschl- Teller potentials for which we ap- 
ply only case 1 of Kovacic's algorithm. Case 1 allows us to obtain the enumerable 
set A„ C A, which includes the classical results obtained by means of supersymmet- 
ric quantum mechanics. Cases 2 and 3 of Kovacic algorithm also can be applied, 
but are not considered here. 

Eckart potential: V(x) — 4coth(x) + 5, x > 0. 

The Schrodinger equation Hty = A 'J is 

&l^= (4coth(ar) +5- A)*. 
By the Hamiltonian change of variable z = z(x) = coth(a;), we obtain 

a(z) = (1-z 2 ) 2 , y(z) = 4z + 5, V(z) 



(z + l)(z- 1)2" 

Thus, K = C(z) and K = C(coth(x)). The algebrized Schrodinger equation Hty 
A$ is 

(1 - z 2 ) 2 <9 2 $ - 2z(l - z 2 )d z ty - (4z + 5 - A)$ = 
and the reduced algebrized Schrodinger equation H$ = A$ is 

4z + 4-A 2-£ |-1 -| 1-| 



~ (z-l) 2 (z + l) 2 (z-1) 2 (z-1) (z + 1) 2 (z + 1) 

We can see that this equation could fall in any case of Kovacic's algorithm. 
Considering A = 0, the conditions {ci,C2,ooi} of case 1 arc satisfied, obtaining 

[Vr\-i = [Vr\i = [Vr\oo = a+ = 0, a ± 1 = a x = 1, af = 2, a± = -1. 

By step 2, the elements of D are and 1. The rational function tu for n = and 
for n — 1 must be 

1 -1 
uj = 1 . 

z+1 z-1 

By step 3, we search for a monic polynomial of degree n satisfying the relation 
(|A.1[) . Starting with n = the only one possibility is Po(z) = 1, which effectively 
satisfy the relation (|A.1[) . while P%(z) — z + ao does not exist. We have obtained 
one solution using Kovacic algorithm: 
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this means that S A„. We can obtain the second solution using the first solution: 



0,2 



+ z - 4 - 4 In (z + 1) z - 4 In (z + 1) 



, *cv. 



0.2 



Furthermore the differential Galois groups and eigcnring s for A = arc 
DGal(L /C(z)) = G„, DG&\(L Q /K) = BGal(L /K) = <G {2} , 



£(H) = Vcct 



' (i-z)2 " i " (i-z) 3 ; ' 



£(#) = Vect 1, 



(coth(x) + l) 2 



-d x 



{l-zf 
coth 2 (a;) + 3 coth(x) +2 



(1 - coth(a;)) 2 (l - coth 2 (x))^ (1 - coth(x)) 3 

Now, for A 7^ 0, the conditions {c2,ooi} of case 1 are satisfied: 



Wr\-l = [aA]i = [V^]oo = a+ = 0, = 1, 



a_i = — § — - a i = a — • 



By step 2 we have the following possibilities for n 6 Z + and for A E A: 



A = 4 



(n+l) 



A+— ) n = a+ -cd-a; =-i + &O i A = 4- -?i 2 -2ti, 



A-+-) n = a x - o+ x - a l = 



A_ 



n = — a_i — oti = 



A = 5 - 4 - n 2 , 
A = 5 - A - n 2 . 



Therefore, we have 



A„ C <U 



— n 2 — 2n : n G Z + ^ U 5 ^ — n 2 : 

(n + 1) 2 II n 2 



(n + l)(z + l)(z-l) ' r ™ (n + l)^(z+l)^(z-l)^ 



Now, for A £ A, the rational function uj is given by: 

^ \ , , _ z(n-l)-n 2 -2n-l ^ _ -2z 2 (n- l)+4z(ra+l) 2 + (ra + l) (n 3 +3n 2 +2n + 2) 

A+— ) 
A-+-) 



nz(n + l) + 2 
(n + l)(z + l)(l-z)' 

z(n-2) — n 2 
n(z + l)(z-l) ' 



(ti+lj^z + lpfz-l) 2 



_ -2z 2 (n-2)+4ra^z + n(n 3 -n + 2) 



A ) CJ = nz(n-l)+2 



■.(»+!)(!-«) > 



7i 3 z 2 (re-l) + 4n 2 z + n 3 -?i 2 +4 



where r„ is the coefficient of the differential equation <9 2 $ = r„$. 
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By step 3, there exists a monic polynomial of degree n satisfying the relation (|A.1[) . 

A 1 r) 2 P J- 9 ( z(n-l)~n 2 -2n-l \ a p , 2(l-n) p — O 

A ++-J O z f n + Z y („+l)(2+l)(z-l) J °z r n + ((„+i) 2 (z+l)(z-lK ™ - U, 

A fl2p | ( nz(n+l)+2 \ q p , nfn+1) p _ n 

A_ 



fl 2 P | o g(rt-2)-» 2 *\ a p 2(2-n) S _ n 

"/ U 2 r " + Z ^(z+ljfz-l) j yzr " + n 2 (z+l)(z-l) r " — U ' 



A \ PP-V I o f nz(n-l)+2 \ o p , n(n-l) p n 

A J C z ^„ + 2 ^ n ( 2+1)(1 _ z) J Cz^n + ( z +l)( z -l)^n ~ U - 

The only one case in which there exists such a polynomial P n of degree n is for 
). The solutions of the equation H$ = A$, with A 7^ 0, are: 

A++_) $„ = P„/„$ , $0 = ^ / n = (2-l)W( Z+ l)5U, 

^ ^ n(l-n) -n(Ti + 3) 

A+— ) $„ = P„/„$ , $0 = 1^ /n= 

A_+_) *„ = P„A,*i, $1 = 3=5; /n = (* + l) 2 ^ 2 (Z - 1) ~" 2 t 3 "- 2 , 



+2 -re z +3re-2 



A—) $„ = P«/„$1, *l = f±I /» = (*+!) 2 " (g - 1) 

In any case \t„ = . The case A_| ) includes the classical results obtained 

by means of supersymmetric quantum mechanics. Thus, replacing z by coth(:r) we 
obtain the eigenstates The eigenring s and differential Galois groups for n > 
and A € A„ satisfy 

~DGdl(L\/ K) C G {2m} , DGal(L A /A') C G {2m} , 

DGal(L A /C(z)) = G m> 

dim c(2) £(H + A) = 2, £(H + A) = £(H + A) = Vect(l). 



Scarf potential: V(x) 



sinh~ x— 3 sinh x 
cosh 2 x ' 



The Schrodingcr equation H^> = E^> is 

„, / sinh 2 x — 3 sinh 2: \ 

a 2 * = 5 E *• 

V cosh a; / 

By the Hamiltonian change of variable z = z(.t) = sinh(x), we obtain 

z 2 -3z 



a(z) = l + z 2 , V(z) 



1 + z 2 



Thus, A" = C(z, Vl + £ 2 ) and AT = C(siim(x), cosh(x)). The reduced algebrized 
Schrodingcr equation H$ = A$ is 

„ 2 ( A2 2 - 12z + A - 1\ 
91 = ' 4(z2 + 1)2 ' *■ A = 3-4A. 
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Applying Kovacic's algorithm for this equation with A = 0, we see that it does 
not satisfy case 1. So we consider only A 7^ 0. By conditions {02, 002} of case 1 we 
have 

[v^-i = [y/r\i = Iv^co = 0, a± = ±±^S 

a + = 5 _ 1 a -. = _ij.i a + = 5 , i r _ _ 1 _ i 

i 4 2' i 4^2 i 4^2' i 4 2' 

By step 2 we have the following possibilities for n £ Z + and for A e A: 
A+++) n = a+ - ati - a+ = > /X + r ~ 4 , A = 4n 2 + 16n + 15, 

A+__) n = a+ - cQ - aj" = aA + I+2 , A = 4n 2 -8n + 3, 
thus obtaining 

A„ C {4n 2 + I6n+ 15 : n € Z+} U {4n 2 - 8n + 3 : n £ Z+} . 
Now, the rational function uj is given by: 

A+++) ^2(z 2 + l)' A +-^ 2(z 2 + l)' 

By step 3, there exists Pq = 1 and a polynomial of degree n > 1 should satisfy 
either of the relation (|A.1[) , 

A "\ 32p i 5z-2 a p Ttz 2 («+4)+3z+n 2 +4n-3 g _ n 

A "1 F)2 p I 2-z a p rtz 2 (7t-2)+3z+rt 2 -2w-3 p _ n 
■<H ) °z*n "r" °z^n {z^+Tp ™ — 

In both cases there exists such a polynomial P„ of degree n > 1. Basis of solutions 
{$i )7 i, $2,n} of the reduced algebrized Schrodingcr equation are: 

A+++) $i,„ = P„/„$i,o, $1,0 = (l + z 2 )fe- arctanz , /„ = 1, 

$2, n = to2,0, $2,0 = M+ai|+y+fa , e -«"*an», = 1. 

A+— ) $!,„ = -Pn/„$1,0, ^,0 = ^™'^, /n = l, 

$2,n = Q„?n$2,0, $2,0 = f±p e^ Ct ™ Z J ^^e" 2 " *dz, Q n = 1. 

In both cases W = -^===, but the classical case (see references |16|, I23j ) is 
A^ ), so that replacing z by sinhec and A by 3 — AE we obtain the eigenstates 

The eigenring s and differential Galois groups are 

£(H - A) = £{H - A) = £(H - A) = Vcct(l), 

DGal(L A /fQ = DGal(£ A /^) = DGal(L A /C(x)) = B. 
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Poschl-Teller potential (revisited): V(r) = cos ^ l~^J^ 2 , x > 0. The 
reduced algebrized Schrodinger equation H<t> = E§ is 

o2 /Az 4 - (A + 3)z 2 + 8\ 

Considering A = and starting with the conditions {c2,ooi} of case 1, we 
obtain 

[v^]o = [Vr\-i = [Vr\i = [v^oo = a+ = 0, a" = 1, 



a + , = at — t, a i = a 



-l ~~ "1 ~ 4' 1 -l ~ "l — 4 



i = -3. a cf = 2 > a o = - 1 - 



By step 2, the elements of D are and 1. For the rational function u, we have the 
following possibilities for n = and for n 

) n = 0, cj 

1 ) n = 0, £J 

A 1 ) n = 1, u; 

A 1 ) n = 1, u; 



1: 






5/4 , 




_ ^1 


z+1 T 




z 


-1/4 




_ ^1 


z+1 




z 


5/4 


z^ H 


_ ^1 


z+1 T 




z 


-1/4 




-1 


z+1 




z 



Following step 3, we search for a monic polynomial of degree n satisfying relation 

(jA.lj) . Starting with n = the only possibility for A+^ ) and A^ | ) is Po(z) = 

1, which does not satisfy relation (jA.lj) in both cases, while Pi(z) = z+ao effectively 

does exist, taking ao = — | for A | ) and ao = § for A | ). This way, we have 

obtained two solutions ($1,0, ^2,0) using Kovacic's algorithm: 



This means that £ A n . Furthermore, 

DGal(L /C(x)) = G [4] , DGal(L /K) = DGal(L /K) = e, 

dime £(H) = 2, dim c £(5) = dim c £(H) = 4. 
Now, for A ^ 0, we see that conditions {c2, 001} of case 1 leads us to 



[v^o = [V^-i = [^1 = [V^oo = 0, 



= 2, 



= -1. 
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A 



++++ 



A 



+++- 



A++-+ 
A++— 
A+-++ 



A+-+- 



A, 



A, 



n = a+ > -atx-ai - <4 



4 



n = a' — a_i - a, — a n 



n = at — a_i — at — at 



2 ' 

/A+T-2 
2 ' 

/A+T-5 
2 ' 

/A+T+l 
2 ' 

/A+T-5 

2 ' 

/A+T+l 
2 ' 

/A+T-2 
2 ' 

_ 7A+T+4 



A = An 2 + 32n + 63, 
A = 4n 2 + 877, + 3, 
A = 4n 2 + 20ti + 24, 

A = 477 2 - 477, 

A = 4n 2 + 20n + 24, 
A = 4n 2 - 4?7, 
A = 4n 2 + 877 + 3, 
A = 4?7 2 - I677 + 15, 



obtaining A„ C A a U A& U A c U A^ U A e , where 

A a = {4?7 2 + 32?7 + 63 : n e Z+} , A 6 = {4r7 2 + 8n + 3 : n e Z + } 
A c = {4t7 2 + 20t7 + 24 : 77 e Z + } , A d = {An 2 - An : n E Z+} , 
A e = {4?7 2 - I677 + 15: ne Z+} . 

Now, the rational function uj is given by, respectively: 



A ++++ ) " = f+4 + f^ + §- A 



a ++ - + ) "=m+^+i a ++- 

A+-++) c = ^ + |4 + |, A+_ H 



A+-+) - = ^ + T^ + i A + ___) u 



5/4 


5/4 


-1 


2 + 1 T 


2-1 ~ 


2 ' 


5/4 


-1/4 ^ 


-1 


z+1 T 


2-1 ^ 


2 ' 


-1/4 


f 5/4 


-1 


2 + 1 


' 2-1 ' 


2 ' 


-1/4 


f " 1/4 




2+1 


r 2-1 
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By step 3, there exists a monic polynomial of degree n satisfying relation (|A.1[) : 



A++++ 
A+++- 
A++-+ 
A++— 
A+-++ 



A+-+- 



A+- 



A+— 



8 2 P 

8 2 P 

8 2 P 
8 2 P 
8 2 P 
8 2 P 



(3z+2)(3z-2) , 



d 2 P 



z(z+l)(z 


-1) 


3z +2 




z(z + l)(z- 


-1) 


2 

6z — 3x — 




z(z+l)(z- 


-1) 


3z-2 




2(a+l)(l- 


,) 


6z 2 +3z- 




z(z + l)(z- 


-1) 


3z+2 




z(z+l){z- 


-1) 


3z 2 -i 




z(z + l)(z- 


-1) 


3z 2 -2 




z(z+l)(l- 


,) 



OzPn 
d z Pn 
d z Pn 
d z Pn 

d z P n 



ti(ti+8) 
n(n+2) 



Pn = 0, 



+ l)(l-z) jr " - U > 



nz(n-\-5)-\-G p 
nz(n— 1) — 12 p 



z(n +5) — 6 p _ 



nz(n— 1)+3 p 
z(z+l)(l-z)- r « 

(rt+2) 



£)2p I 3Z--4 a p I run-t-zj p _ n 



9zP„ 



i(4— n) 



■Pn = 0. 



The polynomial P n of degree n exists for A„ G A& with n even, that is, A n = 
{n G Z : 16n 2 + 16n + 3}, for A ++ _ + ) and A + __ + ). Therefore E — E n = {n G Z 
-4n 2 - 4tt,}. 

The corresponding solutions for A n are 



A 



+++-> 



$1,™ = P2n/n$l,0, $ 



1,0 — ; In 



1, #1,0 



A- 



) 



2,71 



Q2n/n$2,0, $2,0 = 



s^ 5 



/n = 1 #2,0 



These two solutions are equivalent to the same solution of the original Schrodingcr 
equation and correspond to the well known supersymmctric quantum mechanics 
approach to this Poschl- Teller potential, I16L I17j . Furthermore, for all A G A„, 

DGal(L A /C(x)) = G [4] , DGal(L A /A") = DG&l(L x /K) = e, 

dim c £(H - A) = 2, dime £{H - A) = dim c £(H - A) = 4. 

6.4. Searching for Potentials From Parameterized Differential Equa- 
tions. In order to search for new potentials using d z , our main object will be the 
family of differential equations presented by Darboux in |20j , which can be written 
in the form 

(6.10) 8 2 y + P8 z y+{Q - XR)y = 0, P, Q, R G K. 

We recall that some Ricmann differential equations, presented previously, corre- 
spond to this kind. 

When we have a differential equation in the form (|6.10[) . we reduce it to put 
it in the form of the reduced algebrized Schrodinger equation H$ = A<£>, checking 
that Card(A) > 1. Thus, starting with the potential V and arriving to the potential 
V, we obtain the Schrodingcr equation Hty = A^P. This heuristic methodology is 
detailed below. 
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(1) Reduce a differential equation of the form (|6.10[) and put it in the form 
H$ = A<&, checking that Card(A) > 1; to avoid triviality, a should be a 
non-constant function. 

(2) Write W = \d z {\na) and obtain V{z) = a(V - d z W - W 2 ). 

(3) Solve the differential equation (d x z) 2 = a and write z = z(x), V(x) = 
V(z(x)). 

To illustrate this method, we present the following examples. 
Bessel Potentials 

• (From Darboux transformations over V = 0) In the differential equation 

we see that A = —71(71 + 1) and a = z 2 . Applying the method, we obtain 
V = fi and V(z) = fiz 2 + \ with z = z(x) = e ±x . Thus, we have obtained 
the potentials V(x) = V(z(x)) = [ie ±2x + \ (compare with [271 §6.9]). 

• (From Bessel differential equation) The equation 

d 2 y + -d z y + ~ 2 U y = 0, n e \ + Z, 

z z A 2 

is transformed into the reduced equation 
~ (n 2 4z 2 + l\ 

We can see that A = — n 2 , a = z 2 , obtaining V = — z 2 — \, V = — z 4 — 
jz 2 + j and z = z(x) = e ±x . Thus, we have obtained the potential 
V{x) = V(z(x)) = -e ±ix - \e ±2x + \ (compare with HU §6.9]). 

We remark that the previous examples give us potentials related with the Morse 
potential, due to their solutions are given in term of Bessel functions. 

We can apply this method to equations such as Whittaker, Hypergeometric 
and in particular, differential equations involving orthogonal polynomials (compare 
with US §5]). 



Final Remarks 

In this paper we gave, in contemporary terms, a formalization of original ideas 
and intuitions given by G. Darboux, E. Witten and L. E. Gendenshtei'n in the 
context of the Galois theory of linear differential equations. We found the following 
facts. 

• The supcrpotcntial is an algebraic solution of the Riccati equation asso- 
ciated with a potential, defined over a differential field. 

• The Darboux transformation was interpreted as an isogaloisian transfor- 
mation, allowing to obtain isomorphisms between their eigenrings. 
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• We introduced the Hamiltonian algebrization method, which in particular 
allows to apply algorithmic tools such as Kovacic's algorithm to obtain 
the solutions, differential Galois groups and eigenring s of second order 
linear differential equations. We applied successfully this algebrization 
procedure to solve problems in Supersymmetric quantum mechanics. 

• We can construct algebraically solvable and non-trivial algebraically quasi- 
solvable potentials in the following ways. 

(1) Giving a potential where, for A = Ao, the Schrodinger equation is 
integrable. After we put A 7^ Ao, we check that the Schrodinger 
equation is integrable for more than one value of the parameter A. 

(2) Giving a superpotential to obtain the potential, after which we check 
whether the Schrodinger equation is integrable for more than one 
value of the parameter A. 

(3) Using parameterized second order linear differential equations and 
applying an inverse process in the Hamiltonian algebrization method. 
In particular, we can use algebraically solvable and algebraically 
quasi-solvable potentials, special functions with parameters (for ex- 
ample with polynomial solutions). 

This paper is a starting point to analyze quantum theories through Galoisian 
theories. Therefore open questions and future work arise in a natural way: super- 
symmetric quantum mechanics with dimension greater than 2, relationship between 
algebraic and analytic spectrums, etc. 

As a conclusion, as happen in other areas of the field of differential equations, 
in view of the many families of examples studied along this paper, we can conclude 
that the differential Galois theory is a natural framework in which some aspects of 
supersymmetric quantum mechanics may appear more clearly. 
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Kovacic in 1986 (see |36j ) introduced an algorithm to solve the differential 
equation d%( = r(, where r £ C(x), see also [H [24]. 125] [62] . We recall this algo- 
rithm to set the notations. 

Each case in Kovacic's algorithm is related with each one of the algebraic sub- 
groups of SL(2, C) and the associated Riccatti equation 



There are four cases in Kovacic's algorithm. Case 1 concerns the case when 
the Riccati equation has a rational solution (or more); cases 2 and 3 concern cases 
when the Riccati equation admits an algebraic solution; case 4 is the non-intcgrable 
case. 
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Our differential equation is given by <9 2 C = r£ where r = f , s, t € C[x]. We use 
the following notations: 

(1) T' is the set of (finite) poles of r: V = {c G C : t(c) = 0}. 

(2) r = r'u{oo}. 

(3) The order o(r c ) of r at c £ V is the multiplicity of c as a pole of r. 

(4) The order o (7*00) of r at 00 is the order of 00 as a zero of r, 
i-e 0(700) = deg(i) - deg(s). 

Case 1. Let [i/rlc (resp. [a/^Ioo) represent the Laurent series (if defined) of \pr 
at c (resp. infinity). For p g F, we will define e (p) G {+, — } below. Finally, the 
complex numbers aJT, a~, aj,, aQ will be defined in the first step. If the differential 
equation has no poles it only can fall in this case. 

Step 1. For each c e T: 
(co): If o (r c ) = 0, then 

|Vrl c = o, af = 0. 

(d): If o (r c ) = 1, then 

[V^] c = 0, af = l- 

(ca): If o (r c ) = 2, and 

r = 6(.t — c)~ 2 + • • • , then 



r n n ± i±vi + 45 

[vrj c = 0, a c = . 

(c 3 ): If o (r c ) = 2v > 4, and 

r = (o (a? - c) -0 + ... + d (x - c)~ 2 ) 2 + b(x - c)" (t,+1) + • • • , then 
[\/r] c = a (x — c) v + ... + d(x — c) 2 , af = - + v^j . 

(001) : If o (r^) > 2, then 

(002) : If o (t"oo) = 2, and r = • • ■ + 6a; 2 + ■ ■ • , then 



W r \ 00 = °> "oo = o • 



(003): If o (Too) = —2v < 0, and 

r = (ax v + ... + df + bx^ 1 + ■■■ , then 

IVrloo = axV + ••• + d > a7lrf a ^ = \ ~ w 



Step 2. Find D 7^ defined by 
D = < n £ Z + : n = 



cer' J 



If D = 0, then we should start with the case 2. Now, if Card(Z?) > 0, then for each 
n e D we search w eC(i) such that 

a; = e (00) [Vr] ^ + £ (e (c) [v^] c + (x c)" 1 
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Step 3. For each n G D, search for a monic polynomial P n of degree n with 

(A.l) dlP n + 2ud x P n + (d x oj + lo 2 - r)P„ = 0. 

If success is achieved then d = P n e-I w is a solution of the differential equation. 
Else, case 1 cannot hold. 

Case 2. 

Step 1. For each c£f' and for oo, we define the sets E c C Z and C Z as 
follows: 

(ci): If o (r c ) = 1, then E c = {4} 

(c 2 ): If o (r c ) = 2, and r = h b(x - c)~ 2 + • ■ • , then 



{2 + fcVl 



E r = 12 + fcVl + 46 : k = 0, ±2 



(c 3 ): If o (r c ) = v > 2, then S c = {v} 

(001) : If o( roo ) >2, thenS^ ={0,2,4} 

(002) : If ( r oo) = 2, and r = ■ ■ ■ + bx 2 + • • • , then 

Eoo = {2 + Wl + 46 : k = 0, ±2| . 

(00 3 ) : If o (r^) = v < 2, then ^ = {v} 

If some E c = % then switch to next case. 
Step 2. Compute D defined by 



D = { n G Z H 



^ ^eoo - ^e c ^ ,Ve p G S p , p G rj . 



If D = 0, then switch to case 3. If Card(Z?) > 0, then for each n G D 7 compute the 
corresponding rational function 9 defined by 



9 ^ 



2 x — c 

Step 3. For each n G D, search (via linear algebra) for a monic polynomial P n 
of degree n, such that 

(A.2) d*Pn + 38d 2 x P n + {3d x d + W 2 -4r)d x P n + (d x 28 + Wd x 8 + 6 s - 4r6 - 2d x r) P„ = 0. 

If P n does not exist, then case 2 does not hold. If such a polynomial P n is found, 
set = + d x P n /P n and let u be a solution of 

lo 2 + c/)uj + - (d x (j> + 2 - 2r) = 0. 

Then Ci = ^ is a solution of the differential equation. 

Case 3. 

5<ep 1. For each c € V and for 00 we define the sets E c C Z and C Z as 
follows: 

(ci): If o(r c ) = 1, then E c = {12} 
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(02): If o (r c ) = 2, and r = b(x — c) 2 + ■ ■ ■ , then 

E c = 1 6 + ky/l + 46 : fc = 0, ±1, ±2, ±3, ±4, ±5, ±6^ . 
(oo): If o (roo) = u > 2, and r = • • • + bx 2 + ■ ■ ■ , then 

{1 2fc 1 
6 + — 71 + 46 : fc = 0, ±1, ±2, ±3, ±4, ±5, ±6 L m G {4, 6, 12}. 

5£ep 2. For m <E {4, 6, 12} successively, compute the set D defined by 



„ „ to 

D = i n e Z+ : n = — 



If Card(D) > 0, then for each n £ D with its respective m, compute the rational 
function 

_ 12^x-c 

and the (denominator) polynomial defined as S = (x — c). 

cer' 

/Siep 5. For each n € D, with its respective m, compute a monic polynomial 
P n = P of degree n, such that its coefficients can be determined recursively by 

P-i = 0, P m = -P, 

Pi_i = -55 x p - ((m - ») 9 X 5 - 50) P - (m - i) (i + 1) 5 2 rP t+1 , 

where i £ {0, 1 . . . , m — 1, to}. If P does not exist, then the differential equation is 
not intcgrable because it falls in Case 4. 

If P exists, a solution of the differential equation is given by C = e^^, where w is a 
zero of the polynomial of degree m 

i=0 v ' 



Case 4. In none of the above holds, then the equation is non-integrable and its 
Galois group is SL(2). 

Remark A.l ( 4 ). If the differential equation falls only in the case 1 of Kovacic's 
algorithm, then its differential Galois group is given by one of the following groups: 

II: e when the algorithm provides two rational solutions. 

12: G[™1 when the algorithm provides two algebraic solutions Ci, C2 such that 

Cr,C 2 " G C(x) and CrStf" 1 i C(x). 
13: <&{ n } when the algorithm provides only one algebraic solution ( such that 

C" S C(a;) with n minimal. 
I4-' G m when the algorithm provides two non-algebraic solutions. 
15: G a when the algorithm provides one rational solution and the second 

solution is not algebraic. 
16: B when the algorithm only provides one solution £ such that £ and its 

square are not rational functions. 
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